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EDITORIAL 


Chers Bronzés, 


J'espère que vos doigts de pieds n’ont pas eû trop de mal pour revenir dans vos 
chaussures. La position en éventail a du bon, maïs il ne faut pas trop en abuser ! 


La nouvelle la plus importante de cette rentrée est sans aucun doute la décision 
des membres du Bureau de ne pas se représenter lors de la prochaine Assemblée 
Générale. Voir l’article en page 2. 


Des nouvelles en provenance des Etats-Unis : JPC Rom pour le HP-71 est 
maintenant diffusé outre-atlantique. Nous sommes même dans le catalogue 
EduCALC ! 


Des nouvelles en provenance de Paris : des bonnes volontés se sont déjà 
manifestées pour reprendre les programmes de l’ancienne Bibliothèque de 
Genève, entreposée dans une cave parisienne. Continuez à nous écrire si cela 
vous intéresse. 


Plus près de vous : ce JPC fait 32 pages. Nous avons passé tout ce qui nous restait 
comme articles. Il n’y a donc rien de prêt pour le mois prochain. Quelques 
articles pour HP-41 sont en cours de frappe, mais cela ne fait pas un journal 
complet. Alors, à vos plumes et à vos claviers. Rappelez-vous que JPC est votre 
journal. 


A bientôt donc, 


Pierre David (37) 
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A PROPOS DU CLUB 


J'ai à vous faire part d’une nouvelle importante 
concernant le Club. 


Le Bureau, lors de sa dernière réunion, a décidé de 
ne pas se représenter lors de l’Assemblée Générale 
qui aura lieu en janvier, et de passer la main à une 
nouvelle équipe. 


Pourquoi cette décision ? Pourquoi nous arrêter alors 
que le Club semble fonctionner parfaitement et 
donner satisfaction à tous ? 


Nous avons pris cette résolution pour un ensemble de 
raisons que je vais essayer de vous résumer. 


Tout d’abord, il faut savoir que tous les actuels 
membres du Bureau sont possesseurs de HP-71. Nous 
avons adopté cette machine à sa sortie, et l’avons vu 
progresser dans les colonnes du Journal. Aujourd’hui, 
le HP-71 n’est plus fabriqué par Hewlett-Packard. 
Certes, le HP-71 figurera en bonne place pendant 
encore longtemps. Mais, pour nous, ce n’est plus 
pareil. 


A l'inverse, le HP-28 connaît un remarquable succès 
dans JPC. Cette machine n'étant pas au coeur de nos 
préoccupations, nous nous sentons un peu étrangers à 
tout ce qui se fait autour de cette petite merveille. 


Ces constatations viennent à un moment crucial pour 
nous. J’esssplique. Lorsque nous avons repris le Club 
en 1986, la situation n’était pas très brillante. Certains 
d’entre vous s’en souviennent peut-être encore, les 
lettres restaient souvent sans réponse, les journaux 
paraissaient plus qu'épisodiquement, et nous n’étions 
guère plus de 100 membres. Notre objectif a été de 
remettre le Club en état de marche. Aujourd’hui, 
nous considérons que cette tâche est accomplie. 


Il est donc temps de laisser la place aux jeunes ! Le 
Club ne nous appartient pas, nous n’avons pas à le 
monopoliser. Nous avons tous vécu une expérience 
passionnante, mais il faut savoir passer la main. J'ai 
personnellement investi cinq ans dans ce Bureau 
(avec une interruption d’un an), et j’aimerais voir des 
nouvelles têtes. 


Comme je l'ai dit, la situation du Club est florissante. 
Les journaux paraissent à la date prévue, les réunions 
ont lieu à l'heure dite, les finances se portent bien, le 
courrier est traité rapidement. 
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Dans le même ordre d'idées, le Club fonctionnait 
autrefois grâce au matériel personnel de ses 
membres, ce qui posait d'énormes problèmes. Le 
Club a donc acheté dernièrement un ordinateur HP 
Integral PC (sous Unix) pour continuer à fonctionner 
indépendamment des personnes qui donnent déjà 
suffisamment de leur temps. Cet ordinateur facilite 
grandement la fabrication de JPC, la gestion des 
adhérents et autres tâches annexes. C’est un auxiliaire 
précieux et indispensable. Grâce aux outils que nous 
avons développés, il nous fait gagner un temps 
incalculable ! 


Comme vous pouvez le constater à travers ces lignes, 
la santé du Club est bonne. Mais en quoi consiste 
exactement notre activité au Bureau ? Concrètement, 
celle-ci se résume à quelques tâches bien précises. 


La fabrication de JPC, grâce à l’ordinateur du Club, 
est rapide et simple. Rapide, car une journée suffit à 
quelqu'un d’habitué pour mettre en forme tous les 
articles. Simple, car les outils développés sont d’un 
emploi élémentaire et gèrent automatiquement toute 
la présentation de JPC. La fabrication du Journal est 
une activité extrêmement intéressante. 


Le courrier prend, quant à lui, environ 3 heures par 
semaine. C’est une tâche importante, tant par son 
volume que par son intérêt. Le courrier est en 
quelque sorte une vitrine du Club. 


La comptabilité suppose une heure ou deux de travail 
par semaine. Rassurez-vous, notre comptable actuel 
n'a aucune qualification particulière pour cette 
activité. Il suffit d’être méthodique et de ne pas se 
tromper dans ses additions ! 


Il reste encore tous les travaux qui font la vie 
quotidienne du Club. Je cite, pêle-mêle, les visites à la 
Poste, les dépôts de JPC chez notre imprimeur, les 
achats de timbres, les mises sous pli de ZPC, les 
réunions de Bureau, etc. 


Vous avez sans doute compris ce que nous attendons 
de vous. L'Assemblée Générale de janvier élira un 
nouveau Bureau. Nous souhaitons que quatre (au 
moins) parmi vous se manifestent pour le constituer. 


Ceux d’entre vous qui veulent participer devront se 
faire connaître avant le 15 novembre, date limite. La 
transition sera ainsi grandement facilitée. Nous 
aurons le temps de vous expliquer plus en détail le 
fonctionnement du Club. 


Je vous ai résumé le plus succinctement possible les 
raisons qui nous ont amené à ne pas nous 
représenter, ainsi que l’état actuel du Club. C’est 
maintenant à vous de le prendre en charge. Notre 
tâche se termine en janvier. La vôtre commencera 
alors, et vous conduirez le Club comme vous 
l’entendrez. 


Un dernier point : nous restons, bien sûr, fidèles à nos 
chères petites machines. Confier la responsabilité du 
Club à une nouvelle équipe nous laissera plus de 
temps pour mijoter de beaux articles ! 


Quoiqu'il en soit, je vous assure que s'occuper d’une 
association demande du travail, mais apporte 
beaucoup de satisfactions. Durant ces cinq années, j'ai 
rencontré des personnalités vraiment très différentes 
et j'ai beaucoup appris à leur contact. Croyez-moi, la 
vie associative est une expérience passionnante et 
irremplaçable. 


J'espère que vous serez nombreux à répondre à notre 


appel, 


Pour le Bureau, 
Pierre David (37) 





PAIEMENTS HORS DE FRANCE 


Cet article s’adresse à nos membres hors de France 
qui sont amenés à régler leurs cotisations ou toutes 
autres choses par Eurochèques. 


Il s’avère malheureusement que notre banque, ainsi 
que d’autres que nous avons consultées, prélèvent lors 
de l’encaissement de ceux-ci des commissions très 
importantes. A titre d’exemple, nous avons 
récemment reçu un Eurochèque d’un montant de 
75 F sur lequel nous avons touché moins de 20 F. 


Nous appliquons déjà le même tarif à tous nos 
adhérents (même s’ils habitent hors de France), alors 
que les frais d’affranchissement sont différents. 
Actuellement, nos adhérents étrangers coûtent de 
l'argent au Club, ce qui n’est pas très sain. 


Nous avons donc décidé d’ouvrir un Compte Postal 
(numéro 18 823 40 C à Paris) car il nous semble que 
c’est la méthode la moins chère aussi bien pour vous 
que pour nous pour nous payer de l’étranger. 


À partir de ce jour, nous ne pourrons plus accepter de 
paiements que par les moyens suivants : 


-virement sur notre 
18 823 40 C à Paris, 

- chèque en francs français tiré sur une banque 
française. 


Compte Postal numéro 


C’est la seule méthode que nous ayons trouvée pour 
ne pas augmenter votre cotisation. 


Janick Taillandier (246) 





S.O.S. 


Pierre Cazal 

Résidence I Minelli - Bât. H 
Ville di Pietrabugno 

20200 Bastia 


Après avoir voulu introduire dans mon HP-71 un 
programme assez long (15 cartes), le lecteur a refusé 
de lire la huitième carte et répondait obstinément r/u 
Error. Dans pareilles circonstances, existe-t-il un 
moyen de récupérer le contenu des cartes, sans 
retaper le programme comme j’ai dû le faire ? 





COURRIER DU COEUR 


Janick Taillandier 
335 rue Lecourbe 
75015 Paris 


Vend : 
Un module HHP 32 Ko Ram + 32 Ko Eprom : 1000 
FF. 


Pierre Colignon 

24 bis route de Corbeil 
91360 Villemoisson sur Orge 
Tél : 69 04 00 47 
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Vend : SPIE-CAPAGG 
M. Augustin Viguier 


IDS I, II et IT pour le HP-71 : 250 FF. Parc St Christophe Newton 3 
95854 Cergy-Pontoise 
Imprimante ThinkJet HP-2225B : 2000 FF. Tél: (1) 34 22 70 86 
Cherche : 


ee Un HP-97 en bon état. 


Christian Gajac 
Canère à Espiet 
33420 Branne 


Achète : 
HP-71 + module HP-IL + module 
Forth/Assembleur. 


Johan Pyfferoen 
Steenbeukenstraat 3 
8460 Koksijde 
Belgique 


Vend : 
Un MLDL (4 Ko Ram et 4 Ko Rom) de Microbaud 
avec Assembler 3 + manuels + listings de Rom : 
1700 F. 


Claude Becker 
Thierry Laget 

27 avenue Verdi 
59110 La Madeleine 
Tél : 20 55 92 14 


Vendent : 
Un HP-71B, un module HP-IL, un lecteur de cartes 
et une imprimante ThinkJet. 


Eric Gengoux 

8 rue de Furstenberg 

75006 Paris 

Tél : (1) 46 33 65 79 ou (1) 42 85 68 23 


Vend: 

Un HP-28C très peu servi, non modifié : 900 F, une 
imprimante thermique HP-IL HP82162A : 1000 F, un 
lecteur de cassettes HP82161A : 1300 F, une mini 
table traçante Sicape HP-IL : 1200 F, module Ram 4 
Ko pour HP-71 : 200 F. 
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HP28 


P. Heïlbronn Ajustement de courbes d 
S. Lalande Modulation, démodulation 9 
P. Courbis & S. Lalande Solution du grand jeu 10 
P. Courbis Choc en retour 10 
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AJUSTEMENT DE COURBES 


Les lignes qui suivent ont été inspirées par la lecture 
des articles de Simon Bradshaw The HP-27S et 
Curve-fitting on the HP-28 parus dans la revue anglaise 
DataFile et de Wlodek Mier-Jedrzejowicz. 


Simon Bradshaw souligne, pour le HP-27$, la facilité 
d’entrer les données grâce à la structure de liste. Les 
programmes présentés ici transposent ces avantages 
sur la HP-28. On peut ainsi entrer des séries de une 
ou de deux valeurs, afficher des 3x, y et, par touche 
soft, 2xy. On peut interrompre l’entrée des données 
par ICONTJ et la reprendre à tout moment. 
Remarque : dans les lignes qui suivent les touches 
(softkeys) apparaissant dans les libellés de menus 
seront désignées, comme sur les nouvelles HP-22 et 
: HP-32, encadrées par € et }, les signes d’accolade. 


Les programmes présentés ici permettent d’afficher 
en même temps le coefficient de détermination des 
quatre courbes linéaire, exponentielle, logarithmique et 
puissance en indiquant la courbe la mieux ajustée, 
puis de calculer pour chacune, au choix, les 
coefficients a et b, ainsi qu’une estimation de y pour 
tout x. 


1) L’ajustement des quatre courbes est obtenu par 
des appuis de touche similaires. Afin d'éviter toute 
confusion avec la touche tLN], le libellé bilitère du 
programme logarithmique (Ln) comporte une lettre 
minuscule. Pour obtenir une uniformité d’aspect, les 
libellés des trois autres programmes sont bilitères et 
comportent une minuscule. À tout moment, la 
matrice statistique courante DAT peut être éditée. 


2) A la seconde ligne du menu, se trouve la touche 
qui permet de vider de son contenu la matrice 
statistique {ms} et d'indiquer le numéro d’ordre de la 
valeur à entrer dans la matrice (procéder, pour entrer 
deux valeurs, selon la méthode décrite dans le manuel 
d'application pour HP-28C Probability and Statistics 
au chapitre Curve Fitting), puis par appui de la 
touche {5}. Pour plus de commodité, la seconde ligne 
de menu comporte aussi le libellé €F1X). Il est possible 
de fixer de 1 à 11 décimales. Si l’on tape 0 €F1X), l’on 
obtient le mode sTD (appelé ALL sur la 27S). 


Par exemple : taper .123456789123 [ENTER] 
2 {FIX} 0.12 
11 {FIX}  0.12345678912 
O €FIX} .123456789123 


Après l'entrée des données, deux possibilités se 
présentent : 
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a- Vous souhaitez connaître la courbe la mieux 
ajustée. Appuyez sur €R2). Le programme affiche (et 
imprime) les quatre coefficients de détermination et 
indique le meilleur d’entre eux. Si vous voulez obtenir 
des détails sur la courbe, appuyez sur [USER] pour 
aller à la première ligne du menu et choisissez le 
libellé correspondant au meilleur ajustement indiqué 
précédemment; vous pouvez naturellement choisir 
lun des trois autres modèles. Le programme 6 est 
constitué par le programme DELI modifié tel qu'il 
figure dans le manuel HP (page 17) cité ci-dessus. 


b- Vous ne souhaitez pas connaître quelle est la 
courbe la mieux ajustée; appuyez alors sur [USER], 
puis, à la première ligne du menu, choisissez le 
modèle adéquat. Cependant, si vous vous ravisez et 
jugez qu’il vous est utile de savoir quel est le meilleur 
modèle, allez en seconde ligne du menu et appuyez 
sur {R2). 


3) Le libellé Y permet d'accéder par une touche 
partagée entre les quatre modèles (avec les 
indicateurs spécifiques pour chaque modèle) au 
programme de prévision de la valeur de y en fonction 
de la valeur de x (interpolation), dans la mesure où 
les limites suivantes sont respectées : 
- exponentielle : tous les points y > 0 
- logarithmique : tous les points x > 0 
- puissance : tous les points x et y > 0 


La mémoire du HP-28C étant limitée, il est 
recommandé que les programmes ci-dessous décrits 
soient les seuls contenus dans la machine. En outre, il 
est indispensable d’entrer ces programmes dans 
l'ordre indiqué : dans le menu, de bas en haut, de 
droite à gauche. 


MENUS 


1)] ZDAT LI EX LN  PW OO Y | 


2)[ Z MORE ZXY m  R2 FIX | 
|====sssssssssssssssssssssssssssssss | 
3)| A 8 Ç D , E F | 
Josseesesscesesassesmesenssses-am2ze | 
4)  G ZPAR | 


Seules les première et seconde lignes du menu sont 
d’un usage pratique, la troisième et la quatrième 
lignes étant réservées aux routines de service. 


MODE D'EMPLOI 


1) Réinitialiser la mémoire. Désactiver COMMAND, UNDO 
et LAST. 


Routines de service 


2) Entrer >PAR: 
1:23 47 
ZPAR [STO] 


3) Entrer G (routine d’aide à la saisie de données, 
pour les programmes {mR} et {MORE }) : 
« 
CLEAR ‘Zy= ! 1Zx= 1! TOT ARRY+ DROP 
TA 
IF FS? 
THEN 1 - SWAP 1 - 
ELSE SWAP 
END -STR ROT SWAP + 1 DISP 
+STR + 2 DISP 
HALT 


IF FS? 

THEN NZ 1 N> 
START Z- 
NEXT DROP 2 DEPTH ROLL 
START 2+ 
NEXT 

END 

7 CF 8 CF " FIX" 3 DISP 1 WAIT KILL 
» 


G [STO] 


4) Entrer F: 
« 
“ #u NE 1 + 
8 
IF FS? 
THEN 1 + 
END 
DEPTH 3 
IF == 
THEN DUP 1 - 4 ROLLD 
END 
STD +STR " 2" + + 3 DISP 
» 


F CSTO] 


5) Entrer E : 
« CORR LR SWAP » 
E [STO] 


6) Entrer D : 
« SWAP + ROT +STR + DEPTH ROLLD » 
D [STO] 


7) Entrer c: 

Taper 134 [CHR] (niveau 1: "=" 

Utiliser tEDIT) et entrer en ligne de commande : 
« DUP2 # "MMS IFTE » 

C [STO] 


8) Entrer 8 : 
« CORR SQ DEPTH PICK STOZ » 
B [STO] 


9) Entrer A : 


« 


« 
1 RCLZ SIZE LIST+ DROP2 
FOR j 
RCLZ j n 2 +LIST GET LN 
RCLZ j n 2 +LIST ROT PUT STO 
NEXT 
» 
» 


A [STO] 
Deuxième ligne du menu 


10) Entrer fix: 
« DUP FIX 1 IF < THEN STD END » 
fix [STO] 


11) Entrer R2 : 
« 
CLEAR RCLÈ 
B 
2AB 
1AB 
1 À 2 A CORR SQ 
€ 4 1 } +ARRY DUP 
STO> 
ARRY+ DROP 
MAXZ 
C npyit 
C “bn! 
6 "Ex"! 
C "Lin D 
DROP STOZ 
& DISP 3 DISP 2 DISP 1 DISP 


Oo © © 


» 


R2 [STO] 


12) Entrer ns: 
« 
7 SF CLZ CLEAR 
[CT 1 1 1] 'XDAT' STO 
Mt ww 2 DISP: 1 DISP 
" FA 
NZ STD +STR 
" ?" 
+ + 3 DISP 
» 
m [(STO] 


13) Entrer >xy : 


« 


CLEAR RCLÈ DUP 
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TRN 

OVER * 

ARRY+ 

DEPTH ROLL STOZ 

DROP2 ROT DROP2 
» 


Z2xy [STO] 


14) Entrer MORE : 
« F 8 SF G » 
MORE [ESTO] 


15) Entrer 2: 
« F € 2 } +ARRY Z+ G » 
Z [STO] 


Première ligne du menu 


16) Entrer * : 
« 
1 1F FS? 
THEN LN 
END PREDV 
2 IF FS? 
THEN EXP 
END 
» 


Y [STO] 


17) Entrer Pu : 
« 
1 SF 2 SF 
RCLZ 
1A2AE EXP 
& ROLL 
STOZ 
» 


Pw [STO] 


18) Entrer tn: 
« 
1 SF 2 CF 
RCLZ 
1AE 
4 ROLL 
STOZ 
» 


Ln [STO] 


19) Entrer Ex : 
« 
1 CF2SF 
RCLZ 
2 A E EXP 
4 ROLL 
STOZ 
» 


Ex [STO] 
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20) Entrer Li : 
«1CF2CFE » 
Li [STO] 


21) Afin de voir tous les libellés de menu à leur place, 
on pourrait entrer : 1 'ZDAT' [STO], mais le fait de 
taper [s] crée la variable DAT. 


Remarque : Avec une machine dotée de plus de 
mémoire (comme la HP-28$S), on pourrait insérer 
avant la deuxième ligne du menu d’autres modèles à 
étudier. 


Référence : Curve fitting for programmable calculators 
de William M. Kolb. 


EXEMPLE 


(DataFile) : [USER] [NEXTI (m} 
Entrer les données suivantes en deux étapes : tout 
d’abord les trois premières données, puis les deux 


dernières : 


x/ 0.72 1.31 1.95 2.58 3.14 


y/ 2.16 1.61 1.16 0.85 0.50 
Presser ts}. Entrer les trois couples de données : 
.72,2.16€23 1.31,1.61€23 1.95,1.16€2} [CONTI] 


Effectuer des calculs. A présent, vous décidez d’entrer 
les deux couples restants : (MORE) 


2.58,.85€2} 3.14,.5€2)  [CONT] 


Supposons que deux décimales sont désirées pour le 
résultat de R2 : 2 €FIX>. Presser €R2). 


Choisissez, par exemple, l'option a) : 
Pressez €R2) : 
Li 0.99 
Ex 0.98 
Ln =0.99 
Pw 0.90 
Choisissons LN : [USERI€LN} -0.99 (-0.995) (meilleur 
modèle) : b = -1.10, a = 1.85 
Presser [CLEAR] 


Quelles sont les valeurs de y pour x=1.5 etx=2 ? 
1.5 [Y] 1.40 
2 [VIe 1:09 
Choisissons EX : presser €EX} -0.99 (-0.990) b = -0.58, 
a = 3.45 


Presser [CLEAR] 


Quelles sont les valeurs de y pour x=1.5 etx=2 ? 
151" 1.44 

2 CY] 1.08 
Presser [CLEAR] 


Supposons qu’à présent on souhaite connaître les 
quatre coefficients de détermination avec trois 
décimales : 


Presser (NEXT], 3 [FIX]€R2} 
Li 0.986 
Ex 0.980 
Ln =0.989 
Pw 0.900 


EXEMPLE 


(HP Direct Sales Consultant pour HP-18C, page 11). 
Prévisions en utilisant la courbe exponentielle. Vos 
ventes pour un nouveau produit sont les suivantes 
pour les six premiers mois : 


Mois Ventes 
Juin 317 
Juillet 525 
Août 483 
Septembre 566 
Octobre 727 


Novembre 909 


Utilisant le modèle exponentiel, prévoyons les ventes 
pour décembre. Presser {ms}. Entrer : 


317C2} 5252} 483(2} 566€2} 727C2} 9092) 


On pourrait noter y=3527. Presser [CONTI, puis 
2 {FIX} [USER] €Ex} 


Les résultats s’affichent : 


R=0.95 b=0.18 a=293.96 
7 {Y} : 1057.79 Ventes prévues pour décembre. 


Vous souhaitez vérifier quel est le meilleur modèle 
(coefficients de détermination) avec trois décimales. 
[NEXT] 3 CFIX}CR2) 

Li 0.900 

Ex 0.897 

Ln 0.803 

Pw 0.881 


Il aurait été plus conforme à la réalité de choisir le 
modèle linéaire. On aurait alors trouvé pour 


décembre une prévision de ventes d’un montant de 
952.73. 


NB. : Si l’on souhaite imprimer la matrice statistique, 
il est conseillé d'utiliser une HP-28$, en mode sTp, 


avec appui simultané des touches [ON] et EL]. 


Philippe Heïlbronn (233) 





MODULATION DEMODULATION 


Ce programme a un effet saisissant : il incrémente ou 
décrémente les fréquences d’un beep entre deux 
bornes. Il est difficile de décrire le résultat avec des 
mots. Essayez le ! En changeant les paramètres vous 
pouvez faire dévaler les fréquences en une 
bourrasque de sons ou bien faire tomber des 
gouttelettes de cristal sur une plaque de marbre. 


Vous fixez : 

- al : la fréquence de fin en hexadécimal, 

- a2 : la fréquence de début, 

- 43 : la durée de chaque beep, et 

- ad : l’incrémentation. 

Vous remplacez CA et 889 par EA et 88D aux endroits 
marqués d’une étoile pour passer respectivement 
d’une incrémentation à une décrémentation. Il est 
clair que dans un sens il faut que al soit supérieur à 
a2 et dans l’autre, l’inverse. 


Voici le programme et en marge ladaptation pour 
HP-28C (version 1BB) : 


HP-28S HP-28C (1BB) 
69c20 prologue assembleur 
95000 Longueur 
8F18050 sauvegarde de (8F2EE40) 
DO D1BD 
3400010  C-<01000 A1 
06 RSTK=C 
3401000  c=00010 A2 
06 L1 RSTK=C 
D7 D=C 
3405000 Cc=00050 A3 
8F3DA22 exécute Le beep (8FC6101) 
07 C=RSTK 
DA A=C 
3401000  C=00010 A4 
CA A=A+C 
* (EA A=A-C pour 
la décrémentation) 
07 C=RSTK 
D5 B=C 


JPC 57 Page 9 


06 RSTK=C 

Dé C=A 

889 ?B>=C 

*  (8BD ?B<=C pour 

La décrémentation) 

6D GOYES L1 

8F8B050 récupération (8F91F40) 
DO D1BD 

142 A=DATO A 

164 DO=D0+5 

808c PC=(A) 


Vous l’assemblez et l’exécutez directement. 

Je laisse libre cours à votre imagination pour changer 
les paramètres al a2 a3 et ad il n’est pas nécessaire 
d’avoir des talents pour la musique seulement un peu 
d’intuition. 


L’essayer c’est l’adopter. Bon amusement ! 


Sébastien Lalande (442) 


SOLUTION DU GRAND JEU 


Voici le listing désassemblé du programme paru dans 
JPC 54. 


69c20 CON(5) #02C96 prologue de routine LM 
53000 CON(5) #00053 longueur de code 

808F INTOFF inhibe Les interruptions 
133 AD1EX sauve D1 dans A 
1F410F4  D1=(5) #4F014 

300 LC(1) 0 

150 DATI=C 1 écrit 0 en #4F014 
1F4100C D1=(5) #C0014 

15D0 DAT1=C 1 écrit 0 en #C0014 

131 D1=A 

142 A=DATO A fin de routine 

164 DO=D0+ 5 

808C PC=(A) 

Explications 


Le clavier de la HP-28 fonctionne par interruptions : 
INTOFF inhibe donc le clavier. 


Lorsque #4F014 est a zéro (pour la HP-28C, #C0014 


sinon) larrêt système et les auto-tests sont inhibés 
(voir JPC ST). 
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Le but du programme est donc : 


1) d'empêcher l'entrée de caractères au clavier 
2) supprimer larrêt système pour éviter que l’on 
puisse revenir a l’état normal 


Comment s’en sortir 


La seule touche non inhibée est [ON] et on sait 
qu'éteindre puis rallumer la HP-28 remet #4F014 
(#C0014) à une valeur non nulle (voir JPC SI) : il 
faut donc réussir à éteindre la HP-28... 

Pour cela la manière la plus simple est d’attendre son 
auto-extinction. 


Pour débloquer le programme il fallait donc : 
1) attendre l’auto-extinction, 

2) rallumer la HP, et 

3) effectuer un arrêt système (alors possible) 


Simple, non ? 


Paul Courbis (392) 
Sébastien Lalande (442) 


CHOCS EN RETOUR 


A propos de Particle Le facteur est constant (acte I) 
paru dans JPC 56 : 


Philippe Heïlbronn fait les remarques suivantes : 

- on ne peut purger les noms spéciaux de la manière 
habituelle, 

- on ne peut faire de ORDER, et 

- on ne peut pas visiter ces variables. 

En fait, ces trois actions sont toujours possibles : il 
suffit, grâce à ASSEMBLEUR, de placer le ou les noms 
dans la pile, puis d’effectuer les commandes de 
manière classique. 


Le programme ASSEMBLEUR proposé ne fonctionne que 
sur HP-28C, sur HP-28$ remplacer #50000 par #D0000. 


Dans l’article Jeu intégral de Hervé Thévenon (même 
JPC), on peut simplifier le INKEY proposé en le 
remplaçant par : 

« DO UNTIL KEY END » 


Ceci est suffisant... 


Paul Courbis (392) 


HP75 


E. Gengoux 
R. Atlan 


E. Gengoux 
E. Gengoux 


Passerelles entre HP-75 et PC 
Droites de Tukey 


Programme "CSV2VC75" 
Programme "VC752CSV" 


12 
13 


27 
27 
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PASSERELLES ENTRE HP-75 
ET PC SOUS VISICALC 


Plusieurs logiciels permettant à présent d'échanger 
des données sous forme de fichiers texte entre nos 
machines (HP-71 et HP-75) et le monde MS-DOS, je 
vous propose un point rapide sur la question, et des 
utilitaires permettant d’effectuer de tels échanges, 
avec un double objectif : éviter les acrobaties, et aller 
au-delà des simples fichiers texte, mal adaptés à des 
données tabulées. 


Pour cette fois-ci, nous examinerons le cas du HP-75 
et de ses fichiers VisiCalc. Les deux programmes 
proposés, Vc752csv dans le sens HP vers PC et 
csv2vc75 dans l’autre sens, utilisent l’un des formats 
d'échange de fichiers tableurs les plus courants, le 
format C.S.V. (Comma-Separated Values). Ceci 
‘ permet d’envisager des échanges bidirectionnels avec 
la majorité des tableurs disponibles sur PC, mais aussi 
avec d’autres logiciels, statistiques en particulier, qui 
utilisent ce format de données. Que les possesseurs 
de HP-71 ne se désolent pas : dans un prochain 
article, ils trouveront un logiciel analogue, qui leur 
permettra des transferts de fichier DATA (ou, 
beaucoup mieux, de fichier HPAF) vers fichier CSV 
et vice-versa, d'où des possibilités d'échanges de 
données entre HP-71 et HP-75, le plus simplement du 
monde, avec une disquette ! 


LES POSSIBILITES DES PRODUITS 
DU MARCHE 


Voici quelques temps, pour faciliter les échanges 
entre HP-110 ou HP-150 et PC, Hewlett-Packard 
avait introduit la carte d'interface HP82973A, qui se 
monte dans un slot court à l'intérieur de tout 
compatible possédant de tels ports d’extension. Cette 
carte n’est pas chère (environ 1700 F HT), mais son 
interface-utilisateur est limité et pas très clairement 
documenté. Puis est apparu chez HP un logiciel 
beaucoup mieux fait, mais malheureusement limité en 
pratique au seul HP-71, permettant en particulier 
l'échange direct de fichiers entre PC et HP-71 via la 
boucle HP-IL et l’interface précitée, et leur stockage 
sur le disque dur du PC, avec en sous-produit la 
possibilité d'utiliser le PC comme écran-clavier et 
imprimante du HP-71 : logiciel HP-LINK, référence 
HP82477. 


Du côté du monde extérieur, deux produits sont 
apparus quelques mois après ces annonces HP, et 
constituent à notre avis un réel bond en avant du 
point de vue facilité d'utilisation de la carte interface 
HP82973A. Il s’agit, d’une part, du logiciel LINK de 
Southern Software Inc. (à ne pas confondre avec 
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HP-LINK !), d'autre part d’un utilitaire écrit par 
CMIT. pour étendre les possibilités de ses 
remarquables RAM-cassettes, et qui permet à un 
compatible MS-DOS muni de la carte interface HP 
précitée, de lire directement tout support magnétique 
HP (ainsi, une disquette de 75 sur un 9114A ou B 
connecté au PC) et même, dans les cas favorables, de 
convertir au passage les données en format ASCII. 
Cet utilitaire s'appelle CMT-1010 Compatibility 
Software, et je vous engage vivement à l’acquérir chez 
CMT ou chez EduCalc si vous possédez la carte 
HP82973A ! 


Pour résumer, dans les cas favorables c’est à dire 
fichiers ASCII ou fichiers DATA de structure simple, 
un PC pourra lire directement sur disquette les 
données provenant d’un HP-71 ou d’un HP-75. En 
sens inverse, il faudra utiliser l’un des deux logiciels 
HP-LINK ou LINK puisque ni le HP-71, ni le HP-75, 
ne peuvent relire de support formatté sous MS-DOS, 
et on sera limité aux seuls fichiers de type texte dont il 
ne faudra pas oublier de renuméroter les lignes, dans 
le cas du HP-75 ! 


LES PROGRAMMES DE CONVERSION 


Le problème posé est donc d'utiliser un type de 
fichier texte particulier, qui conserve la structure du 
fichier de données dont il est issu. Or, justement, les 
fichiers C.S.V., pour lesquels chaque enregistrement 
est transformé en une ligne de texte (champs séparés 
par des virgules, et zones alphanumériques encadrées 
par deux guillemets), répondent convenablement à 
cette exigence, avec les deux avantages 
supplémentaires de pouvoir être édités du moins, tant 
que la longueur de ligne reste raisonnable, et de 
n'occuper en mémoire qu’une place réduite 
(beaucoup plus, en tout cas, qu’un fichier DATA, 
SDATA, DIF...). 


En sens inverse, la difficulté réside dans la 
reconstitution des différents champs en utilisant les 
virgules sauf celles éventuellement présentes dans les 
champs alpha eux-mêmes ; c’est la raison d’être des 
guillemets ! 


Avec le VisiCalc du HP-75, dont la puissance est tout 
à fait remarquable, on évite au passage la difficulté 
liée aux fichiers DATA (HP-71 et HP-75), pour 
lesquels la détection des fins de ligne ou 
d'enregistrement n’est, en accès séquentiel, pas aisée 
lorsqu'on ne connaît pas a priori sa structure, surtout 
lorsqu'elle n’est pas fixe ! Dans ces fichiers, point de 
MAXCOL ni d’accès direct au niveau du champ avec 
GETLABELS$ et GETVALUES, hélas ! 


Je vous laisse découvrir les deux programmes : ils 
peuvent traiter n'importe quel fichier VisiCalc, avec 
ou sans label utilisateur, et quelle qu’en soit la taille. 
Bien entendu, les formules ne sont pas transférées : 
seules les valeurs et les étiquettes le sont. 


L’isolement commence à être rompu, non ? A 
suivre... 


Post-scriptum : pour être complets, notamment dans 
le cas où, possédant un portable genre Toshiba, 
Zénith ou Sharp, vous ne pouvez pas utiliser 
facilement une carte d’extension du type de la 
HP82973A, il existe un autre logiciel de 
communication, avec des possibilités comparables à 
celles de LINK, mais qui ne fait appel qu’au port 
série : 71B-TALK (de Interactive Concepts, chez 
EduCalc). Mais, dans ce cas, il n’est évidemment plus 
question de lire directement des disquettes HP sur le 
PC : par malheur, le procédé d’enregistrement 
physique utilisé par HP dans le 9114 n’est pas le 
même que celui utilisé par les autres fabricants de 
lecteurs de disquettes 3} pouces. Dommage !.. 


Eric Gengoux (108) 





AJUSTEMENT LINEAIRE 
DROITE DE TUKEY 


Partons d’une distribution statistique (x, y;), 
c’est-à-dire d’un tableau à deux colonnes et à 7 lignes. 
Sur un graphique, cette distribution prend la forme 
d’un nuage de points. Notre objectif est de construire 
un "résumé" de cette distribution sous la forme d’une 
droite qui passe au plus près des points du nuage. Il 
nous faut donc déterminer la pente a et l’ordonnée à 
l'origine b de la droite souhaitée. 


La méthode d’ajustement la plus connue est la 
régression linéaire fondée sur la méthode des 
moindres carrés. L'objectif est de minimiser la 
somme des carrés des écarts entre les points du nuage 
et la droite de régression. La régression linéaire ne 
requiert que des calculs algébriques simples et trouve 
naturellement sa place dans le schéma théorique 
gaussien. 


Malheureusement, la régression linéaire n’a pas une 
qualité précieuse en statistique : la robustesse. Les 
coordonnées de tous les points interviennent dans la 
détermination de la pente et de l’ordonnée à l’origine 
de la droite de régression, donc un point déviant pèse 
lourdement sur les paramètres de la droite. La 


s 


"déviance" peut être due à une fluctuation aléatoire 
(bruit ou erreur de mesure par exemple) ou être 
significative (individu étranger à la population) ; dans 
les deux cas, le statisticien souhaite que cette déviance 
ne pèse pas sur les caractéristiques décrivant la 
population étudié, mais que les points déviants soient 
isolés, que les écarts qu’ils présentent soient ensuite 
analysés de façon fine. 


Nous connaissons tous une caractéristique peu 
robuste : la moyenne arithmétique d’une distribution. 
Cette moyenne varie dès qu’une des valeurs de la 
distribution change. Au contraire, la médiane est une 
caractéristique robuste. Rappelons que la médiane 
s'obtient en prenant "la valeur du milieu" de la 
distribution rangée dans l’ordre croissant. Si le 
nombre de points est pair, nous prenons comme 
médiane la moyenne arithmétique des deux valeurs 
du milieu. 


Pour revenir à notre problème d’ajustement linéaire, 
nous allons apprendre à construire la droite dite de 
TUKEY, du nom du statisticien américain qui la 
inventée. Cette construction repose sur la médiane et 
nous fournit, bien sûr, la pente et l’ordonnée à 
l’origine de la droite de TUKEY. Dans cet article, 
nous essaierons d’en exposer simplement la méthode. 


La droite de TUKEY et la méthode des trois groupes 


Quand nous voulons étudier comment y varie avec x, 
nous construisons un diagramme y en fonction de x. Si 
le nuage de points est allongé selon une droite, 
celle-ci représente un bon résumé de la relation de y 
avec x. Si le nuage de points est incurvé, transformer y 
ou x ou les deux variables simultanément (par 
exemple, effectuer une transformation logarithmique 
ou prendre l’inverse) peut être une bonne façon de le 
redresser, et nous sommes alors ramenés au cas 
précédent. Nous supposerons ici que le nuage ne 
présente pas de courbure particulière. 


Comme il est probable que quelques points du nuage 
seront déviants par rapport à la tendance générale, 
notre façon de travailler visera à ne donner qu’une 
importance minime à ces points. Ainsi ces points 
déviants généreront des résidus (écarts entre eux et la 
droite de TUKEY) importants. L'étude des résidus 
nous permettra peut-être alors de trouver des 
explications plus fines les concernant. L'idée de base 
est de diviser l’ensemble des x points qui composent 
la distribution statistique en trois groupes égaux, de 
caractériser chaque groupe à l’aide de médianes et de 
construire la droite de TUKEY à partir des points 
médians. 
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Construire les trois groupes 


Nous commençons par trier les valeurs x; de sorte 
qu’elles soient rangées par ordre croissant. A partir 
de ce tri, nous divisons les 7 points en trois groupes, 
de tailles aussi voisines que possible : celui de gauche, 
celui du milieu, celui de droite. 

Quand il n’y a pas d’ex-aequo dans les x;, le nombre 
de points dans chaque groupe dépend du reste de la 
division de n par 3. Habituellement, nous répartissons 
les données entre les trois groupes comme suit : 


Groupe n = 3k n=3%k+1 n=3%k+2 
Gauche k k k+1 
Milieu k k +1 k 
Droite k k k+1 


Points médians 


Pour chacun des groupes ainsi construits, nous 
déterminons les coordonnées de son point médian en 
déterminant d’abord la médiane des x, puis, 
séparément, la médiane des y;. Nous appelons Xg Et Vos 
les coordonnées du point médian du groupe de 
gauche, x,, et y, celles du point médian du groupe du 
milieu et enfin x, et y, celles du point médian du 
groupe de droite. La figure 1 montre les points 
expérimentaux, le découpage en groupes suivant les x; 
et les points médians dans un exemple portant sur 9 
points. Comme illustre la figure, il est possible 
qu'aucun des points médians ne soit un point 
expérimental car les médianes suivant les x et celles 
suivant les y sont déterminées séparément. Dans bien 
des cas, nous pourrons observer néanmoins que les 
points médians coïncident avec des points 
expérimentaux. L'usage de la médiane garantit ainsi 
la robustesse de la droite de TUKEY construite à 
partir des points médians, pour autant que le nombre 
de points dans chaque groupe ne soit pas trop petit. 


Pente et ordonnée à l’origine 


Il nous faut maintenant déterminer la pente a et 
l’ordonnée à l’origine b de la droite de TUKEY : 


ÿ = a*x + b 

l'accent circonflexe nous rappelle que cette équation 
nous permet d’estimer la valeur prise par y pour une 
valeur donnée de x. La pente a est donnée par la 


formule : 


a = (ya * Ya) / Kg * Xg) 
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En utilisant (7 Y9) et (Kw Ya) ainsi, nous répondons à 
deux contraintes : calculer la pente de la droite de 
TUKEY sur un intervalle assez grand, garantir la 
robustesse si nous avons assez de points dans les 
groupes de gauche et de droite. L’ordonnée à 
l’origine est donnée par la moyenne : 


b= COYg ” a*xg) + Cm * XX) + Yq * 2*X4)) / 3 


L'usage des médianes garantit encore ici la robustesse 
du résultat. 


Exemple 


Dans une étude en 1953, Greenberg a mesuré l’âge et 
la taille de 18 enfants d’une école privée urbaine. 
Nous avons reproduit ses données dans le tableau 
ci-dessous. 


Source : B.G. Greenberg, The use of analysis of 
covariance and balancing in analytical studies, 
American Journal of Public Health, 43, 692-699. 


Enfant Age Taille 
(mois) (cm) 
1 109 137.6 
2 113 147.8 
3 115 136.8 
4 116 140.7 
5 119 132.7 
6 120 145.4 
7 121 135.0 
8 124 133.0 
9 126 148.5 
10 129 148.3 
11 130 147.5 
12 133 148.8 
13 134 133.2 
14 135 146.7 
15 137 152.0 
16 139 150.6 
17 141 165.3 
18 142 149.9 
Nous pouvons construire le nuage de point 


représentant les 18 enfants étudiés (cf. figure 2). Ce 
nuage ne présente aucune courbure prononcée, aussi 
nous ajustons une droite aux données pour expliquer 
comment la taille varie en fonction de l’âge dans ce 
groupe d’enfants. Les enfants 13 et 17 semblent être 
déviants par rapport à la tendance générale et nous 
les étudierons de près au fur et à mesure de nos 
travaux. 


Comme 18 est divisible par 3 et qu'il n’y a pas 
d’ex-aequo en x, nous construisons 3 groupes de 6 
points chacun. Les points médians sont alors : 


Xg 5 Yg) = (115,5 ; 139,15) 
CXm ? Ym) = (127,5 ; 147,9 ) 
CXq ï Yg) = (138  ; 150,25) 


Vous noterez que le nombre des points de chaque 
groupe étant pair, nous avons à prendre la moyenne 
du troisième et du quatrième de chaque groupe pour 
déterminer les médianes. 


La pente de la droite de TUKEY est alors : 
a =0,4933 et son ordonnée à l'origine est 
b = 83,118 


La méthode des moindres carrés appliquée à la 
distribution ci-dessus donne : a = 0,511 et b = 79,7 


Les valeurs obtenues par les deux méthodes sont donc 
voisines, ce qui est le cas lorsque la distribution 
étudiée ne présente pas de points déviants. Il faut 
souligner à ce propos que, en plus de son manque de 
robustesse, la régression linéaire classique présente 
un autre défaut souvent négligé : les formules qui 
donnent la pente a et l’ordonnée à l’origine b ne sont 
valides que si les points (x, y;) sont connus 
exactement ; dans le cas le plus fréquent de résultats 
expérimentaux, il faut prendre en compte l'existence 
d'erreurs et utiliser des méthodes plus complexes 
(voir bibliographie). 


Comment programmer la méthode de TUKEY ? 


L'affaire se décompose en modules suivant la 
description que nous venons d’en faire. 


a) module de saisie où il faut déterminer la taille n du 
tableau de rangement des x;, y; où l’on saisit et valide 
les valeurs. 


b) module de découpage en 3 groupes où nous avons 
trois cas possibles : 

- n divisible par 3 

- reste de la division par 3 égal à +1 

- reste de la division par 3 égal à + 2 


c) module de calcul des points médians, où il faut 
appliquer systématiquement un sous-programme de 


calcul de médiane 


d) module de calcul de la pente et de l’ordonnée à 
l'origine et d’affichage des résultats. 


Essayez-vous y ! 


Source de ce document 


J.D. Emerson et D.C.Hoaglin, Resistant lines for y 
versus x extrait de D.C. Hoaglin, F. Mosteller et J.W. 


Tukey, Understanding Robust and Exploratory Data 
Analysis, J.Wiley, New-York, 1983. 
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and Computing of exploratory Data Analysis, Duxbury 
Press, Boston, 1981. 

F. Mosteller et J.W. Tukey, Data analysis and 
regression (a second course in statistics), 


Addison-Wesley, Reading, 1977. 


-sur la régression linéaire classique tenant compte 
des erreurs expérimentales 


KS. Krane et L. Schecter, Regression line analysis, 
American Journal of Physics, 50, 82, (1982) 


J. Orear, Least squares when both variables have 


uncertainties, Am. J. of Phys., 50, 912 (1982). 


Robert Atlan (285) 
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fig.1 - Points expérimentaux (x) et points médians (0) dans un 
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fig.2 - Taille en fonction de l'âge de 18 enfants. 
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fig.3 - Résidus après ajustement d'une droite de TUKEY. 
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AH ! VOUS ECRIVEZ 


Vous vous sentez en verve, mais vous ne savez pas sous quelle forme "l’équipe de 
rédaction" souhaite recevoir votre prose. C’est ici que se trouvent les réponses à 
vos questions. 


Dans la mesure du possible, vous devez nous envoyer vos écrits sur support 
magnétique (carte, cassette ou disquette). Soyez sans crainte, nous vous 
retournerons vos biens après copie. 


Si vous ne pouvez pas utiliser de support magnétique, ou ne pouvez vous rendre 
aux réunions, alors et alors seulement faites le sur papier. 


Que ce soit sur une feuille de papier, ou sur support magnétique, ne dépassez 
pas 50 caractères par ligne. 


Pour nous épargner du travail, insérez dans votre texte les commandes de 
formattage suivantes (et non les commandes du formatteur HP) : 


Au 


centre un titre, par exemple : 
TITRE 


"\"(CHR$(92)) marque le début et la fin d’un paragraphe. Par exemple : 


\Début de paragraphe exprimant le contenu de vos idées qui, même si vous en 
doutez, intéressera certains des membres du Club. Surtout si vous vous sentez 
débutant. Les articles pour débutants écrits par des débutants sont ceux qui 
manquent le plus. Fin de paragraphe.\ 


N'oubliez pas de mettre les accents. Utilisez le jeu de caractères Romans. Les 


possesseurs de HP71 utiliseront les redéfinitions de touches ci-dessous, ainsi que 
le fichier CHARLEX listé dans le coin des Lhex. 


Jean-Jacques Dhénin (177) 


DEF KEY 'fW!, CHR$(197); (é) 
DEF KEY ‘fE', CHR$(193); (é) 
DEF KEY 'fR', CHR$(201); (è) 
DEF KEY 'fY!, CHR$(203); (ù) 
DEF KEY ‘fU', CHR$(195); çà) 
DEF KEY ‘fl', CHR$(209); ci) 
DEF KEY 'f0', CHR$(194); (ô) 
DEF KEY ‘f/!, CHR$(92); (\ 
DEF KEY 'fA!, CHR$(192); (à) 
DEF KEY 'fS', CHR$(200); (à) 
DEF KEY ‘fD', CHR$(205); (8) 
DEF KEY ‘fJ!', CHR$(207); (ü) 
DEF KEY 'fK!, CHR$(221); «i) 
DEF KEY ‘f*1, CHR$(124); €p) 


DEF KEY 'fC!, CHR$(181); (ç) 


PPC PARIS SE REUNIT 
UNE FOIS PAR MOIS 


Comme vous le savez peut être déjà, PPC Paris se réunit une fois par mois, en 
plein coeur de Paris. Amenez votre matériel, votre bonne volonté et vos idées ! 
Plus vous en apporterez, et plus vous en trouverez chez vos collègues de PPC. 


Ces réunions se déroulent de manière très libre, aucun ordre du jour, discussion 
ou autre n'étant imposé, Un membre du bureau est toujours présent. Ainsi, si 
vous désirez remettre votre article tout frais au Journal, si vous avez des 
suggestions à faire, si vous voulez vous procurer des anciens numéros de JPC, ce 
sera en principe toujours possible. 


Si donc cela vous intéresse, n'hésitez plus un seul instant, venez nous rejoindre 
tous les premiers samedis de chaque mois (sauf en période de vacances 
scolaires) au : 


Centre de Jeunesse et de Loisirs Jean Verdier 
11 rue de Lancry 
75010 Paris 


et en montant au deuxième étage, vous entendrez des éclats de rire et des 
discussions passionnées vers la salle 215. Attention, toutefois, de venir entre 16 et 
19h. 


Pour l’accès en métro, trois possibilités s’offrent à vous : 
- Métro Strasbourg Saint Denis : 

Sortie porte St Martin / Bd St Denis, coté pairs 

- Métro République : 

Sortie Bd St Martin, coté pairs 

- Métro Jacques Bonsergent : 

Sortie Bd Magenta, coté impairs. 


Ah, j'oubliais ! JPC est (souvent) distribué en avant première lors de ces 
réunions... À bon entendeur, salut ! 


Les dates des prochaines réunions sont : 
Samedi 17 septembre 1988 

Samedi 1°’ octobre 1988 

Samedi 5 novembre 1988 

Samedi 3 décembre 1988 

Samedi 7 janvier 1989 

Samedi 4 février 1989 

Samedi 4 mars 1989 

Samedi 6 mai 1989 

Samedi 3 juin 1989 


Pierre David (37) 


NOUS EN AVONS 


La coopérative du Club vous propose : 


- de lecteurs de cartes magnétiques pour HP-71, neufs, dans leur boîte d’origine, avec 5 cartes magnétiques, 
pour 500 F (port compris), 

- des anciens numéros de JPC, au prix de 40 F + 7,40 F de frais d’affranchissement, 

- d’une année complète de numéros de JPC (février à janvier) pour 300 F (offre spéciale) port compris, 

- des I.D.S. du module Forth / Assembleur (listing interne commenté par HP) pour 250 F (port compris), 

- des VASM pour HP-41 (listings des Roms internes commenté par HP) pour 300 F (port compris), 

- de manuels de service du HP-41 au prix de 75 F (port compris), 

- de manuels de service du HP-75 au prix de 75 F (port compris). 


En outre, le module JPC Rom pour HP-71 est disponible. Vous nous adressez votre Eprom CMT (de 
préférence 64 Ko), et nous la programmons suivant une des options ci-dessous : 

- JPC Rom + Manuel, pour 600 F, 

- JPC Rom + Manuel + vos propres programmes, pour 800 F. 


Si vous souhaitez des renseignements complémentaires, n’hésitez pas à nous contacter. 


VOUS EN VOULEZ 
Nom : 
Prénom : 
No de membre : 
Adresse : 
Commande : 

Qté Prix Unitaire Prix Total 
le TS Re D M mt 
anciens numéros de JPC x 47,40 FF 
année complète de JPC x 300 FF 
1.D.S. du module Forth x 250 FF 
V.A.S.M. x 300 FF 
Manuel de service pour HP-41 x 75 FF 
Manuel de service pour HP-75 x 75 FF 
JPC Rom + Manuel x 600 FF 
JPC Rom + Manuel + vos propres programmes x 800 FF 
Actualisation Eprom x 150 FF 
A A Ce à ni de nl nf Bla EU ILE 


Préciser éventuellement Les 
numéros de JPC commandés : 





PPC PARIS 


Association régie par la loi de 1901, enregistrée 
à Paris le 2 décembre 1982 sous le numéro 82/3240 


BULLETIN D’ADHESION 
ON EX AS DEL A ES AN 3 4 D PA NOR PF PT A A NP 1 
Prénom |_[_1_|_1_1_1_1_1_1_1_1_1_1_1_1_1LI Date de naissance |_|_1/|_1_1/1_1_| 
A PP PP DEEE PB PE 1 ei 4 A A HE UN A GA EN 
LS SL ER AE Lt IE Ad Pl 
LISTE ERA LES OR ES TE RES RUE ET ASP ENLA 
LEE ET PORTER ET EAN EE ER IEEE TA 
Code Postal |_[_[_[_1_1 Ville [_1_1_1_121IUITUITLILICI ILLICITE 
| PET RCI EE) 


Profession 
Intérêts 


Matériel HP en votre possession 


Autre matériel informatique 


Comment avez-vous connu PPC Paris ? 


Que recherchez-vous au sein de PPC Paris ? 


La Loi No 78-17 du 6 janvier 1978 relative à l'informatique, aux fichiers et 
aux libertés, garantit à toute personne justifiant de son identité un droit 
d'accès et de rectification auprès des services où organismes chargés de mettre 
en oeuvre des traitements informatiques comportant des informations nominatives 
Le concernant. 


Je souhaite adhérer au club PPC Paris conformément aux statuts de 
l'Association. Au mieux de ma connaissance, je déclare avoir Le droit de 
fournir tous Les programmes et informations que je vous enverrai (sans 
enfreindre des obligations de secret à l'égard d'autres personnes ou 
organismes) pour publication dans Le Journal de liaison, sans obligations ni 
responsabilité d'aucune sorte (en cas d'utilisation frauduleuse) de La part des 
dirigeants de PPC-Paris. 


Date |_|_|/1_1_1/191_ 11 
Signature, précédée de La mention "Lu et Approuvé!! 
Le montant de La cotisation s'élève à 350.00 F pour un an. 


Etudiants: 300.00 F (justificatif indispensable) 


Paiement à l'ordre de "PPC Paris! par chèque bancaire ou virement postal (CCP 
No 18 823 40 C à Paris). Ne pas utiliser d'Eurochèques. 


Veuillez envoyer toute correspondance à : 
PPC Paris, BP 604, 75028 Paris Cedex 01, France 


ASSEMBLEUR 


J. Dardennes 
L. Istria 
L. Guillou 


BASIC 


P. Aspero 
L. Istria 
E. Gengoux 


E. Gengoux 
E. Gengoux 
P. Aspero 
L. Istria 

E. Gengoux 


LE COIN DES LHEX 


Interpolation linéaire 
Evaluation polynômiale (acte IT) 
Affaire à suivre. à la trace 


Le facteur est premier 
Tri sur fichier texte 
Passerelles HP-71 et PC 


Programme "CSV2VC75" pour HP-75 
Programme "VC752CSV" pour HP-75 
Programme "PF" 

Programme "TRI" 

Programme "TXTPC71" 


18 


21 


23 
23 
26 


27 
27 
28 
29 
29 


31 
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INTERPOLATION LINEAIRE 


C'est avec enthousiasme que je vous demande de 
renouveler mon abonnement à PPC. 


J'en profite pour vous envoyer le source d’un petit 


C=R2 
GOSBVL DV2-15 
GOSUB  SVR2R3 
GOSUB RESTD1 
GOSUB  AD32D1 
D1=D1+ 16 
GOTO FIN 
NIBHEX 8888855 


Lex sans prétentions, mais qui me rend de grands  LIiNe GOSUB COMMON 
services, car je suis élève en deuxième année de D1=D1+ 16 
l'ENSAM. CD'1EX 
GOSUB DIFF 
Il s’agit de deux fonctions relatives aux interpolations C=R3 4 
linéaires. Elles s’intitulent LIN et INVLIN. D=C 4 
C=R2 4 
Utilisation : GOSBVL MP2-15 
GOSUB  SVR2R3 
Yo = LIN OX4, Y4s X2s Y2r Xo) GOSUB  RESTD1 
GOSUB  AD32D1 
Xo = INVLIN CX4, Yqs X2, Y2, Yo) GOSUB  AD32D1 
GOTO FIN 
J'espère que ce petit Lex pourra rendre service à  DIFF  D1=C 
d’autres, même si on ne peut le comparer à des GOSBVL POPIR 
merveilles telles STRUC2.... C=D 
D1=C 
Longue vie à PPC C=A 4 
GOSBVL POPIR 
Jean Dardennes (369) GOSBVL SPLTAC 
P=C 15; 
LEX ‘LIN! 9%P= O0 
ID #5E GOYES NEG 
MSG 0 POS P= 0 
POLL 0 GOTO Ok 
POPIR EQU  #0E8FD NEG  P= 9 
SPLITA EQU  #OC6BF ok CP 15 
SPLTAC EQU  #0C934 GOSBVL AD2-15 
AD2-15 EQU #0C363 RTN 
MP2-15 EQU #OC43A SVR2R3 R2=A W 
DV2-15 EQU  #OC4AC C=B M 
URES12 EQU  #0C994 R3=C M 
FUNCD1 EQU  #2F8C0 RTN 
FNRTN1 EQU  #0F216 RESTD1 D1=(5) FUNCD1 
ENTRY  INVLINe A=DATI À 
CHAR _#F D1=A 
ENTRY LINe RTN 
CHAR #F AD32D1 D1=D1+ 16 
KEY  'INVLIN! D1=D1+ 16 
TOKEN 3 RTN 
KEY ‘LIN! SUB32D1 D1=D1- 16 
TOKEN 2 D1=D1- 16 
ENDTXT RTN 
NIBHEX 8888855 COMMON GOSUB  PUSHD1 
INVLINe GOSUB COMMON GOSUB  RESTD1 
GOSUB AD32D1 D1=D1+ 16 
CD1EX CD1EX 
GOSUB DIFF DEC  W 
C=R3 W CD1EX 
D=C GOSUB  AD32D1 
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CD1EX 
GOSUB DIFF 
GOSUB  SVR2R3 
GOSUB RESTD1 
GOSUB AD32D1 
CD1EX 
D=C W 
CD1EX 
GOSUB 
CD1EX 
GOSUB DIFF 
C=R3 W 
D=C W 
C=R2 W 
GOSBVL DV2-15 
GOSUB  SVR2R3 
GOSUB  RESTD1 
CD1EX 
D=C W 
CD1EX 
GOSUB 
RTN 
PUSHD1 CD1EX 
D1=(5) FUNCD1 
DATI=C A 
CD1EX 
RTN 
FIN GOSBVL POPIR 
GOSBVL SPLITA 


AD32D1 


AD32D1 


CER3 
D=C W 
CER2 WW 


GOSBVL AD2-15 
GOSBVL uRES12 
GOSUB RESTD1 
GOSBVL FNRTN1 
END 





EVALUATION POLYNOMIALE 
(ACTE Il) 


Dans le numéro 50 de JPC, page 28, Jack Helay nous 
fournissait une fonction d'évaluation de polynômes. 
Pour intéressante qu’elle était, cette fonction péchait 
quand même par un coté. Que faire si, dans un 
programme, il était impossible de connaître à priori, 
le degré du polynôme ? 


Afin de remédier à ce problème, il était nécessaire 
d'écrire une fonction admettant, non pas un nombre 
variable de paramètres, mais un tableau. 


Chose est faite. La syntaxe de cette nouvelle fonction 
est : HORNER(T,X), À étant la variable, et T étant le 
tableau dans lequel sont stockés les coefficients du 
polynôme. 


Pierre David m'avait conseillé de faire en sorte de 
jumeler les deux fonctions. On aurait ainsi pu écrire 
les coefficients du polynôme, soit sous forme de 
tableau, soit in extenso, comme dans la fonction de 
Jack. 


Hélas, trois fois hélas, lors de l'écriture d’une 
fonction, on doit faire précéder le point d’entrée d’un 
certain nombre de quartets qui indiquent au 71 quels 
sont les paramètres valides. le deuxième bit de chacun 
de ces quartets indique la possibilité d’avoir un 
tableau. S'il est à un, l’entrée d’un tableau est valide. 
Dans le cas contraire, elle ne l’est évidemment pas. 


Ce bit ne pouvant pas prendre 2 valeurs en même 
temps (ce serait idéal non ?) il est impossible d’avoir 


le choix entre les deux formes. 


Les heureux possesseurs de HP-71 ont donc à leur 
disposition deux fonctions : INTERP et HORNER. 


Pour le moment, HORNER ne traite que les variables 
réelles (REAL) et réelles courtes (SHORT). On ne peut 
encore effectuer des calculs sur les nombres 
complexes. Je n’ai pas traité les variables entières 
(INTEGER). 
A bientôt. 

Laurent Istria (3) 


LEX "HORNER 


=id EQU #5D 


FNRTN1 EQU #0F216 
ARGERR  EQU #0BF19 
POPIR  EQU #OE8FD 
SPLITA EQU #OC6BF 
STABT  EQU #0D3D9 
MPY EQU #OECBB 
MP2-15 EQU #OC43A 
AD2-15 EQU #0C363 
RCCD1  EQU #0D3F5 
URES12 EQU #0C994 
FUNCD1 EQU #2F8C0 
CON(2) =id 
CON(2) #A 
CON(2) #A 
CON(5) 0 
NIBHEX F 
REL(4) 1+TxTbSt 
CON(4) 0 
CON(5) 0 
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CON(3) 
REL(5) 
CON(1) 


TxTbSt 


TxEn01 


CON(1) 
NIBASC 
CON(2) 
NIBHEX 


CTxEn01)-(TxTbSt) 


Horner 
#F 


11 
1HORNER! 
#A 

1FF 


ARR HE HR ee ee A ee he he eh HR AR Re 


* 


* 


HORNER 


* But: Implémente La fonction HORNER 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


Horner 


Æ # * * *#  * 


ok 


Entrée: 
- P 
- D1 


0 


Stack pointer 


- D0O= Program counter 


Sortie: 
- D1 mis à 
- sort par 


jour 
FNRTN1 


Appel: POPIR, SPLITA, STAB1, MPY, MP2-15 
AD2-15, RCCD1, URES12, FNRTNI 


ARGERR 


Abime: A,B,C,D,P,RO,R1,R3, FUNCD1 


Niveau : 3 

Historique: 
88/06/28 : 
88/06/29 : 
88/07/01 : 


NIBHEX 
GOSBVL 
GOSBVL 
GOSBVL 
SETHEX 
D1=D1+ 
* 
A=0 
A=A#1 
P= 
LC(1) 
B=C 
C=DATI 


Voici un dope 


[addrt |dim2|di 
5 4 4 


?B=C 
GOYES ok 
B=B+1 
?2B=C 
GOYES 
GOVLNG A 
ST=C 


L.I. Conception et codage 


L.I1. Documentation hor20 
L.I1. Traite Les Short hor10 
KARRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRARURRARMMUÉ AR 
8A22 1 tableau et un nombre 
POPIR Le nombre 
SPLITA format 15 chiffres 
STAB1 dans RO R1 
16 D1 * Le dope vecteur du 
tableau 
A préparation au problème 
A de OPTION BASE 
0 
#B type short hor40 
P 
W 
vecteur 
mlo|diT| 
111 
hor60 
P 
P type real boucle 
P 
ok 
RGERR autre type 
Type du tableau dans ST 


CSRC 
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A=A-C 
B=A 
P=C 
D1=D1+ 
ADIEX 
D1=(5) 
DAT1=A 
A=0 
CSR 
A=C 
CSRC 
CSRC 
CSRC 
CSRC 
CSR 
R3=C 
* 
D=C 
C=0 
C=D 
?P= 
GOYES 
A=0 
A=A+1 
GONC 
A=A+B 
C=C+B 
GOSBVL 
CR3EX 
CSRC 
CSRC 
CSRC 
CSRC 
CSRC 
D1=C 
D1=D1- 
?2ST=1 
GOYES 
D1=D1- 
GOSUB 
CR3EX 
C=C-1 
C=C-1 
GONC 
A=DAT1 
GOSBVL 
GOTO 
CR3EX 
SETDEC 
GOSBVL 
GOSBVL 
GOSBVL 
C=B 
D=C 
C=A 
GOSUB 


FUNCD1 
A 


W 
A 
A 


hor40 
read 


W 

W 
horé0 
W 
SPLITA 
Fin 


SPLITA 
RCCD1 
MP2-15 
W 

W 

W 

read 


C(0) OPTION BASE 


1 81:07; Osei 

dans B 

dimension du tableau 
pour le retour 


D1 à sa place 
précaution 

C(A)= dimension 1 
dans A 


pour conserver 
l'adresse du tableau 


précaution de même 
CCA)= dimension 2 
la dimension est 2 ? 


sinon Le première 
vaut 1 
B.EsT: 
rajout 
OPTION 
nombre 


de 1 si 
BASE = 0 
d'éléments 
on récupère Le 
dope vecteur 

et on cherche 
l'adresse du 
tableau 


adresse dans D1 
pour le GOSUB 


pour Le GOSUB 
premier élément 
nombre d'éléments 
pour Le 

compteur 

plus d'un élément 
sinon constante 


évite Les ennuis 
forme 15 chiffres 
X dans C, D 
multiplication 
copie du résultat 
dans C, D 


lecture de l'élément 


GOSBVL SPLITA 
GOSBVL AD2-15 addition 
CR3EX compteur d'éléments 
SETHEX en hexadécimal 
C=C-1 W 
SETDEC pour les autres 
x opérations 
CR3EX remise à sa place 
GONC boucle si il y a encore des 
x éléments on continue 
Fin D1=(5) FUNCD1 on récupère D1 
C=DAT1 A 
D1=C 
XM=0 
GOSBVL URES12 forme 12 chiffres 
GOVLNG FNRTNT et fin 
read ?ST=1 0 short ? 
GOYES short 
D1=D1+ 16 élément suivant 
A=DATI VW lecture de l'élément 
RTN 
short  D1=D1+ 9 élément suivant 
P= 8 pour le Lire 
A=0 w précaution 
A=DAT1 WP lecture de l'élément 
ASRC l'élément n'est codé 
ASRC que sur 9 quartets 
ASRC il faut donc le 
ASRC mettre à La bonne 
ASRC place 
ASRC 
RTN 
END 
AFFAIRE A SUIVRE 
… À LA TRACE 


Voici une fonction Trace. La Trace est la somme des 
éléments diagonaux d’une matrice carrée (première 
dimension = deuxième dimension). 


Les tableaux normaux complexes et réels sont 
acceptés ; mais les entiers, réels courts et complexes 
courts ne marchent pas. 


J'ai repris une partie du programme DEMOCOMP 
de Jean-Jacques Moreau (voir JPC 31, page 20) pour 
le formattage et la sortie des nombres complexes sur 
la math stack. 


Mon programme n’est qu’une ébauche : il peut être 
optimisé et on peut le rendre compatible pour les 
tableaux réduits et entiers. 


Lionel Guillou (326) 

LEX ILEXTRACE! 

ID #5C 

MSG 0 

POLL O0 

ENTRY DEBUT 

CHAR #F 

KEY TRACE! 

TOKEN 101 

ENDTXT 
FNRTN1 EQU #0F216 
URES12 EQU #0C994  Conv. 15-forme/12- forme 
AD2-15 EQU #0C363 Add. 2 nbre 15-formes, 


* respectivement en A,B et C,D;résultat en A,B 


SPLITA EQU #OC6BF  Conv 12-forme(A)/15-forme 
MEMCHK  EQU #012C7 
BSERR EQU #0939A 
EXPR EQU #0OF23C Fin d'une fonction 
SPLITC EQU #0C940  Conv 12-forme(C)/15-forme 
RDATTY EQU #17CCé Message d'erreur 
ERREUR GOTO  POPERR 
NIBHEX A11 Argument :un tableau 
DEBUT 
CD1EX 
R1=C R1 := D1 initial 
CD1EX 
C=0 4 
A=0 W 
A=DATI B Lecture:type de variable 
P= 0 
LCHEX E 
?2A=C P Complexe ? 
GOYES COMP 
LCHEX C 
2A#C P Tableau réel ? 
GOYES ERREUR Ni L'un,ni l'autre:erreur 
GOSUB TEST 


R2=C  R2=nbre de registre à sauter avant 
* de lire Le registre correspondant à l'élément 
* diagonal suivant 
C=DATI À 
* La math stack 
CD1EX 
GOSUB SOMME 
GOSBVL URES12 
GOVLNG FNRTN1 
* pour Le tableau réel 
COMP GOSUB TEST 
C=C+C A R2=nbre de registres à 
R2=C sauter pour Lire Le prochain elt 
* diagonal(rem:il est 2 fois plus important que 
* pour le tableau réel car un nbre complexe est 


Adresse du tableau sur 


fin de la routine 


tableau complexe 


* composé de 2 nbre:réel et imaginaire) 
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C=DAT1 A 
* La math stack 


Adresse du tableau sur 


CD1EX 
GOSUB SOMME 
GOSBVL URES12 Conversion en 12-forme 
SETHEX 
RO=C Stocke somme réel dans RO 
C=R1 R1 = D1 initial 
CD1EX POPERR 
D1=D1+ 11 SOMME 
C=DAT1 A Rappel de l'adresse du 
CD1EX tableau sur La math stack 
D1=D1+ 16 
C=R2 
CSRB divise par 2 
C=C-1 A 
C=C-1 A 
R3=C Nbre d'itérations pour 
* somme 
GOSUB SOMME BOUCLE 
GOSBVL URES12  Conv. en 12 flt 
SETHEX 
D=C W Stocke Le nbre imaginaire 
C=R1 dans R1 
CDEX W et D1 initial dans A 
R1=C 
C=D A 
A=C A 
C=0 A 
P= 0 FIN2 
LCHEX 12 Longueur prise,sur La 
* pile,par La partie im. et l'entête d'un nb comp. 
GOSBVL MEMCHK INCR 
CDEX A Sauvegarde l'éventuel 
* numéro d'erreur ds D(A) 
D1=C Restaure D1 
GONC  STOCKE Y-a-t-il suffisamment de 
* place? 
C=D A Non;récupère Le code 
* d'erreur 
GOVLNG BSERR FIN1 
STOCKE A=RO Partie réel.-> A(W) 
DATI=A W -> Stack; 
D1=D1- 16 
A=R1 Partie im. -> A(W) 
DATI=A W -> Stack; 
D1=D1- 2 
LCHEX O0E en tête de nb complexes 
DATI=C B -> Stack; 
XM=0 Désarme XM 
GOVLNG EXPR Fin de La fonction 
TEST D1=D1+ 2 
A=DATI 5 1ère dimension 
D1=D1+ 5 
C=DAT1 4 2ème dimension 
ASRC 
2A#C A Dimensions égales 


GOYES  POPERR 
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A=0 
ASLC 
C=C-A 
R3=C 
D1=D1+ 
C=C+1 
C=C+1 
RTN 
GOVLNG 
A=DAT1 
B=0 
SETDEC 
GOSBVL 
SETHEX 
C=R3 
C=Cc-1 
GOC 
R3=C 
GOSUB 
C=DAT1 
D=0 
SETDEC 
GOSBVL 
GOSBVL 
GOSUB 
C=R3 
C=Cc-1 
R3=C 
GONC 
C=R1 
CD1EX 
RTN 
C=R2 
SETHEX 
CSL 
D=C 
CD1EX 
C=C+D 
CD1EX 
RTN 
SETHEX 
SB=0 
XM=0 
GOTO 
END 


RDATTY 
y 
W 


SPLITA 


FINI 


INCR 
W 


SPLITC 
AD2-15 
INCR 


BOUCLE 


FIN2 


A:indice La plus bas:0 ou 1 
R3=dim-indice Le plus bas 


1ère lecture 
Conv 15- forme 
R3=compteur de La boucle 


Décrémentation du compteur 


incrémentation de D1 
Lecture du nbre 


Conv 15- forme 
addition:(A,B)=(A,B)+(C,D) 
incrémentation de D1 


décrémentation compteur 


restaure D1 initial 


C*ié(taille d'un registre) 


incrémentation de D1 


LE FACTEUR EST PREMIER 


Je vous propose un petit programme, nommé PF : 
quelque chose comme facteurs premiers, mais en 
anglais. 


Mode d’emploi 


1- Exécutez le programme, 


2- Apparaît N,n= 

Introduire deux nombres séparés par une virgule : 

- le nombre N désignera la valeur à laquelle le HP-71 
devra commencer ses recherches, et le nombre x le 
nombre de nombres premiers cherchés. 


- si le nombre n est égal à 0, l’indicateur binaire 0 
reste allumé, et le HP-71 recherchera tous les 
nombres premiers à partir de la valeur N jusqu’à 
épuisement des batteries. Il appartient à lutilisateur 
de ce programme de fixer l’ordre DELAY à sa 
convenance. 


3- Appuyez sur [ENDLINE] 

Ceux qui disposent d’une imprimante pourront 
aisément rajouter l’ordre nécessaire pour imprimer 
les nombres premiers qu’ils voudront. 

Exemples 

Exemple 1 : 

N,n= 1000,0 

Le HP-71 affiche successivement : 

1009 


1013 
1019 


Exemple 2 : 
N,n= 1000,10 


Après un court moment, le HP-71 affiche les 10 
nombres premiers suivant le nombre 1000. 


Pour les revoir, il suffit de faire RUN 330. 


Philippe Aspero (424) 


TRI SUR FICHIER TEXTE 


Avant d’entrer dans le vif du sujet, permettez moi de 
faire quelques rappels au sujet des tris. 


Tout d’abord, qu’est ce qu’un tri ? C’est une méthode 
qui permet suivant un certain nombre de critères 
d’ordonner les éléments d’une liste donnée. Il est bien 
entendu qu'il existe quasiment une infinité de ces 
méthodes. La meilleur étant évidement le plus rapide, 
le moins gourmande en place mémoire et celle qui 
respecte le plus l’état du tableau tel qu'il était à 
l'entrée, cela pour autoriser le tri suivant plusieurs 
clés. 


Je me bornerai à vous narrer ici la grande aventure de 
4 types de tri. Le tri bulle, le tri par extraction, le tri par 
insertion et le tri Shell. Le programme que je vous 
propose est d’ailleurs bâti autour de cette dernière 
méthode. 


Avant toute discussion, posons déjà une convention. 
Dans tout le reste de cet article, j’appellerai clef, un 
élément précis du tableau que nous trions. clef(i) 
étant bien entendu le "€ élément du tableau. 


Tri bulle 


Commençons par la moins performante de toutes ces 
méthodes, le tri bulle. Pourquoi bulle ? et bien, car 
comme vous l’allez voir, les éléments remontent de 
place en place dans le tableau; on a ainsi l'impression 
de voir se déplacer une bulle d’air qui remonte dans 
l'eau. 


La méthode consiste à comparer les clefs deux à deux 
à partir du bas du tableau, et à échanger les éléments 
s’il y a lieu. Voici un exemple : 


46 46 46 46 46 46 46 * 
1e; 12 2, 22 4 12 1* 12 
53 53 53 53 53 1* 12 12 
15 15 15 15 14:93 53. 53 
33 33 33 DS MS: 15 15 
1 1 1% 1932.93 93: 33 33 
10 8* 8 8 8 8 8 8 
8* 10 10 10 10 10 10 10 
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On remarque que l’élément dont la clef est la plus 
petite remonte de place en place jusqu’en haut du 
tableau. En recommençant cette opération, le 
deuxième élément dans l’ordre des clefs trouvera sa 
place, et ainsi de suite. En effectuant n-1 fois cette 
opération, sur le tableau de 7 éléments, on aura 
ordonné celui-là. 


On peut écrire cette procédure de la façon suivante : 


begin 
for i = haut, bas-1 
for i = bas, haut+1, -1 
if clef(i) < clef(i-1) 
swap(i,i-1) 


end if 
next 
next 
end 
Deux remarques permettent d'améliorer 


sensiblement l'algorithme. Tout d’abord, après j-1 
opérations d'examen du tableau, les j-1 plus petites 
clefs sont triées. Il suffit donc d'examiner les éléments 
dont l’indice est supérieur ou égal à j pour placer la 
jième clef. 


Ensuite, il faut remarquer que plusieurs clefs peuvent 
être placées lors d’une seule opération. Le tableau est 
ainsi susceptible d’être trié avant le (n-1)\è"€ examen. 
On peut donc penser cesser les boucles lorsque le 
tableau est ordonné, et plus spécifiquement lorsqu'il 
est ordonné dès le début des opérations. On va pour 
cela se servir d’un drapeau qui permettra après un 
examen complet du tableau de déceler si celui-ci est 
ordonné. 


Ce drapeau est activé par l'échange de deux éléments; 
on le baisse au début de l’examen, et on le lève lors 
d’un échange. Avant un nouvel examen du tableau, on 
procède à un test, et ne continue que si le drapeau est 
levé. 


On écrit donc la procédure ainsi : 


begin 
j + haut 
repeat 
flag + 0 
for i = bas, j+1, -1 
if clef(i) < clef(i-1) 


swap(i,i-1) 
flag + 1 
endif 
next 
EE 
until j > bas-1 or flag == 0 


end 
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Le nombre maximum d'échanges est de n(n-1)/2, ce 
qui peut être très long pour des tableaux importants. 
On réservera donc ce tri aux tableaux de petite taille. 


Tri par extraction 


Le deuxième type de tri que nous allons voir est le tri 
par extraction. Il consiste à « extraire » l’élément 
ayant la plus petite clef, et à l’échanger avec l'élément 
de gauche du tableau (pour un tableau ordonné de 
gauche à droite). 


Pour trouver le plus petit élément, on explore le 
tableau de gauche à droite, en gardant en mémoire la 
plus petite clef rencontrée. La procédure est la 
suivante : 


begin 
for i = g, d-1 
min + clef(i) 
ind + i 
for j = i+1, d 
if clef(j) < min 
min + clef(j) 


ind + j 
endi f 
next 
swap (i,ind) 
next 
end 


Le nombre d’échanges dans ce cas est, comme pour le 
tri bulle de n(n-1)/2. Mais cette méthode a le 
désavantage de ne pas s’apercevoir lorsqu'un tableau 
est déjà trié. De la même façon, on utilisera ce tri que 
pour les petits tableaux. 


Tri par insertion 


Le troisième type de tri est le tri par insertion. Cette 
méthode consiste à insérer dans une partie déjà 
ordonnée du tableau un élément situé à la limite de 
cette partie, puis à recommencer avec les autres 


éléments : 


gs. éléments ordonnés ----- > 

&r. 12: 225 3771 (397 56! 79. 62 

Î Î 
gauche droite 


Pour insérer l'élément droite dans le sous-tableau, il 
suffit de le sauvegarder en mémoire, de décaler vers 
la droite tous les éléments dont la clef est supérieure 
à la sienne, puis de l’insérer à sa place : 


& Me: 23.37 39 56 79 
<rse clefs. <42.r22-> 


sauvegarde=42 
& 12 23 37 39 62 56 7 


Voici la procédure qui permet de réaliser cette 
opération : 


begin 
for j = g+1, d 

i + j-1 

s + clef(j) 

while i > g+1 and clef(i) > s 
clef(i+1) + clef(i) 
+ 1-1 

end while 


if clef(g) > s 

i + g-1 

clef(g+1) + clef(g) 
endif 


clef(i+1) + s 
next 
end 


J'ai intentionnellement séparé la partie if endif du 
reste du programme. En effet, la condition i z g and 
clef(i) > s de la boucle while peut conduire à une 
erreur. Lorsque i vaut g-1, i > g est faux, mais clef(i) 
n’a pas d’existence. On s’affranchit de ce fait en 
remplaçant i >g par iz2g+1 et en traitant 
particulièrement le cas où à = g. 


On calculera que cette méthode utilise au plus 
n(n-1)/2 déplacements (pour un tableau en ordre 
décroissant), et au moins #-1 (pour un tableau trié). 
Elle est en plus stable, et est donc à préférer pour le 
tri de tableaux de petite dimension. 


Tri Shell 


Voici enfin le tri Shell, du nom de son inventeur. Il est 
de loin le plus rapide des quatre tris que je vous aurai 
proposé ici. 


Shell s’est aperçu que le tri par insertion était rapide 
dans deux cas précis : sur les tableaux déjà ordonnés, 
et sur les tableaux de petites dimensions. D’où son 
idée de diviser les grands tableaux en des sous-unités 
plus petites, et de faire effectuer aux éléments des 
sauts importants à l’aide d’un pas. 


Dans l'exemple suivant, le 2 est placé en 5 tests dans 
le cas d’insertion simple, et en 3 seulement si l’on fait 
effectuer des sauts de deux places aux éléments (pas 
de deux). 


Chacun des sous-ensembles est donc ordonné par la 
méthode du tri par insertion. Dans un tableau 
ordonné aléatoirement au départ, les places d’indices 
(1, 2, 3. pas) contiennent les minima de ces 
sous-ensembles. En faisant varier le pas d’une valeur 
fixée à l’avance jusqu'à 1, on obtient un tableau 
totalement ordonné, mais on travail sur des 
ensembles petits d’une part, et de plus en plus 
ordonnés d’autre part. 

Les valeurs du pas doivent être simples à calculer. 
Peterson et Russel montrent qu’il est très efficace de 
prendre pour P les valeurs p,,, P,n1 Px P, telles 
que : 


P1=1 

Pix1 = Pi t1 

P est le premier élément de la suite tel que 
Pm 2 IP(1/9) 


On arrive alors à cette procédure : 


begin 
p+î 
while p < ip((d-g+1)/9) 
p + 3*p+1 
end while 


repeat 
for k=g,p 
for i = k+p, d, p 
s + clef(i) 
16e #5 
while j 2 k+p and clef(j) > s 
clef(j+p) + clef(j) 
J+jr=p 
end while 
if clef(k) > s 
jeke.p 
clef(ktp) + clef(k) 
end if 
clef(j+p) <- s 
next 


next 
p + ip(p/3) 
until p == 
end 


Délaissant le tri de tableaux numériques, j'ai préféré 
trouver une méthode de tri des fichiers texte. Le 
problème était un peu plus difficile, car la lecture 
d’une fiche ne se fait pas directement, mais par 
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l'intermédiaire d’une variable de stockage, et d’un 
numéro de canal. Les structures du type while et if ne 
pouvaient donc pas être utilisées de façon habituelle. 


J'ai donc été obligé de recourir aux boucles LOOP... END 
Loop, afin de rendre adéquates les structures 
mentionnées plus haut. Je suis alors arrivé au petit 
programme dont vous trouverez le listing plus loin, à 
sa place habituelle. 


Il est bien entendu que la simplicité du programme 
n’a d’égale que sa lenteur d’exécution. Mais n’oubliez 
pas que les fonctions INSERT# et READ#, sont lentes à 
opérer. J'ai entendu parler d’un tri rapide en 
assembleur. Mais peut être n’étaient-ce que des voix ? 


Une fonction simple permet de définir le temps 
moyen de fonctionnement du programme pour un 
fichier contenant n enregistrements. Cette fonction 
est : 

t= n162/11 


Notez que l’on a pris là le cas le plus défavorable, 
puisqu'il s’agit de temps pour des fichiers ordonnés à 
l'envers au départ. On s'aperçoit quand même du fait 
qu’il faut de la patience. 


Et maintenant sortez vos fichiers textes et triez-les 
(Merci M. Hofstadter). 


Laurent Istria (3) 


PASSERELLES HP-71 / PC 


Il y a, dans les logiciels ÆP-LINK, LINK et 
71B-TALK dont nous avons évoqué dernièrement les 
caractéristiques, deux possibilités intéressantes 
vis-à-vis des fichiers texte (ceux qui intéressent 
l’auteur d’articles que je suis !), à savoir le transfert 
de fichiers vers autre chose (WordStar, par exemple...) 
et l'utilisation des ressources du PC (écran et 
imprimante). Malheureusement, certains caractères 
bien utiles en français ne sont pas codés de la même 
manière dans les deux machines, ce qui a des 
conséquences inesthétiques à l'impression... 


Or, dans la dernière version de JPC Rom, il y a une 
fonction absolument géniale (si, si) qui permet 
virtuellement toute transcodification et, luxe suprême, 
travaille directement au niveau du fichier, donc sans 
besoin de dérouler de boucle pour lire une ligne, la 
traiter, la remettre là où il faut, et ainsi de suite... J’ai 
nommé MAP, de Tapani Tarvainen (voir JPC 46). 
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Rappelons-en la syntaxe : 


MAP fichier , chaîne , chaîne2 

MAP fichier , chaînel , chaîne2 , début 

MAP fichier , chaînel , chaîne2 , début , fin 
MAP # canal , chaînel , chaîne2 

MAP # canal , chaînel , chaîne2 , début 
MAP # canal , chaîne , chaîne2 , début , fin 

Le programme proposé peut travailler dans les deux 
sens : de HP vers PC, ou l'inverse. Dans ce dernier 
cas, il permet aussi de donner une équivalence 
approximative pour les fameux caractères 
semi-graphiques. Il crée et utilise deux tables de 
transcodification, une dans chaque sens, et pourrait 
certainement être réécrit sous forme d’un Lex : avis 
aux amateurs ! Nous nous en servons de temps à 
autre au Club, notamment quand un article nous 
arrive sur disquette 5% pouces. Attention, le 
programme ne traite pas les Escape Sequences (ce 
serait possible avec REPLACES, mais dangereux dans 
certains cas...). Il est très, très rapide. 


Eric Gengoux (108) 


Programme "CSV2VC75" pour HP-75 (transfert HP-75 / PC) 


0010 ! CSV2VC75 Utilitaire Xfert 71/PC > 75 
NB: Zones alpha obligatoirement entre "1, 
Numéroter Les Lignes et éliminer Les en-têtes HPAF ou autres 
“descriptor blocks! si présents. 
(c) Eric GENGOUX 13\04\1988 V1.0 
0040 DIM R$(256] 
0050 INPUT “Source CSV !";: F1$ à F1$=UPRC$(F1$) 
0060 IF INCAT(F1$,'T')#2 THEN BEEP à GOTO 50 
0070 INPUT "Dest worksheet !'; W$ @ W$=UPRC$(WS) 
0080 IF INCAT(W$,"W")#1 THEN BEEP @ GOTO 70 
0090 DISP "User headers in CSV [i/6]" 
0100 K$=UPRC$(WKEY$) à IF NOT POS('YN',K$) THEN BEEP @ GOTO 100 
0110 DISP "OK. Wait..." @ IF K$='Y!' THEN K=0 ELSE K=1 
0120 ASSIGN # 1 TO F1$ @ WORKSHEET W$ 9 R=K-1 
0130 ON ERROR GOTO 290 ! Fin de fichier 
0140 READ # 1 ; R$ à R=R+1 à C=K-1 ! Boucle ext.(R) 
0150 C=C+1 ! Boucle int.(C) 
0160 IF R$[1,1]1='#! THEN GOTO 170 ELSE GOTO 230 
0170 P=POS(R$[2],'"1)+1 ! Champ type label 
0180 C$=R$S[I2,P-1] à R$=RS [P] 
0190 IF LENCR$)=1 THEN R$='! ELSE R$=R$(3] ! détection fin Ligne 
0200 IF C+R=0 THEN GOTO 270 
0210 IF (C=0 EXOR R=0) AND (C$=" # OR LEN(C$)=0) THEN GOTO 270 
0220 PUTLABEL C,R,C$ à GOTO 270 
0230 P=POS(R$,',') ! Champ type valeur 
0240 IF P=0 THEN C$=R$ @ R$='! @ GOTO 260 
0250 IF P#0 THEN C$=R$[1,P-1] à R$=R$S[P+1] 
0260 IF C+R#0 AND C$#!' THEN PUTVALUE C,R,C$ 
0270 IF LENCR$)#0 THEN GOTO 150 ! Fin boucle int. 
0280 GOTO 140 ! Fin boucle ext. 
0290 PUTWIDTH 0,7 @ OFF ERROR à ASSIGN # 1 TO !! ! EOF 
0300 BEEP 400,.5 à BEEP 300,.5 à DISP "CSV > VC75 done!" à END 


HONOR HA ARAORORHHOHHHAAAA e A A ACROROR RO A R RRR 


Programme "VC752CSV" pour HP-75 (transfert HP-75 / PC) 


0010 ! VC752CSV Utilitaire Xfert VC75>CSV (71 ou PC) 
0020 ! Ne convertit pas Les formules 

0030 ! (c) Eric GENGOUX 1310411988 V1.0 

0040 DIM R$[256] 

0050 W$=ACTIVES à IF NOT LENCW$) THEN W$='"XXXXXXXX!! 

0060 INPUT "Source US ",W$; W$ @ W$=UPRCS(UW$) 

0070 IF INCAT(W$,"W")#2 THEN BEEP à GOTO 60 

0080 INPUT ‘Destination ","CSVfile!":; F2$ @ F2$=UPRC$(F2$) 
0090 IF INCAT(F2$,"T")#1 THEN BEEP à GOTO 80 

0100 DISP "Row/col headers üeep/ärop!" 

0110 KS=UPRC$(WKEY$) à IF NOT POS("KD",K$) THEN BEEP @ GOTO 110 
0120 DISP "Wait..." 

0130 ASSIGN # 2 TO F2$,TEXT 

0140 IF K$='K' THEN K=0 ELSE K=1 

0150 R$=!'! à WORKSHEET W$ 
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0160 FOR R=K TO MAXROW 

0170 FOR C=K TO MAXCOL 

0180 C$=GETLABEL$(C,R) à T=0 

0189 IF LEN(C$)=0 AND (C=0 OR R=0) THEN C$=! ! à T=0 
0190 IF LENÇ(C$)=0 AND C+R#0 THEN C$=GETVALUES(C,R,0) à T=1 
0200 IF T=1 THEN R$=R$&!,'&C$ @ GOTO 220 

0210 IF T=0 THEN R$=R$&!,"'&C$&'n! 

0220 NEXT C 

0230 R$=R$[21 à PRINT # 2 ; R$ à R$='! 

0240 NEXT R 

0250 ASSIGN # 2 TO '' à BEEP 400,.5 à BEEP 300,.5 
0260 DISP “VC > CSV done! à END 


HR RH he he ee He He He he eee OA A A A AAA ee 


Programme "PF" (calcul de facteurs premiers) 


10 DESTROY ALL @ OPTION BASE 1 à SFLAG O0 

20 INPUT 'N,n=';N,J @ IF J THEN CFLAG O0 à DIM P(J) 
30 IF N<38 THEN 360 

40 M=MOD(N-37,30) @ L=N-M @ IF M<5 THEN 120 

50 IF M<7 THEN L=L+6 à GOTO 130 

60 IF M<11 THEN L=L+10 à GOTO 140 

70 IF M<13 THEN L=L+12 @ GOTO 150 

80 IF M<17 THEN L=L+16 Q GOTO 160 

90 IF M<23 THEN L=L+22 à GOTO 170 

100 IF M<25 THEN L=L+24 Q GOTO 180 

110 L=L+30 à GOTO 190 

120 L=L+4 à GOSUB 200 à L=L+2 

130 GOSUB 200 à L=L+4 

140 GOSUB 200 à L=L+2 
150 GOSUB 200 à L=L+4 
160 GOSUB 200 à L=L+6 
170 GOSUB 200 à L=L+2 
180 GOSUB 200 à L=L+6 

190 GOSUB 200 à GOTO 120 

200 IF NOT MOD(L,7) THEN RETURN 

210 K=11 à M=SOR(L) 

220 IF NOT MOD(L,K) THEN RETURN 

230 IF NOT MOD(L,K+2) THEN RETURN 

240 IF NOT MOD(L,K+6) THEN RETURN 

250 IF NOT MOD(L,K+8) THEN RETURN 

260 IF NOT MOD(L,K+12) THEN RETURN 

270 IF NOT MOD(L,K+18) THEN RETURN 

280 IF NOT MOD(L,K+20) THEN RETURN 

290 IF NOT MOD(L,K+26) THEN RETURN 

300 K=K+30 @ IF K<=M THEN 220 

310 IF FLAG(O) THEN DISP L @ RETURN 
320 I=1+1 @ P(I)=L @ IF I<J THEN RETURN 
330 1=0 
340 1=1+1 @ DISP PCI) à IF I<J THEN 340 
350 END 
360 IF N<3 THEN L=2 @ GOSUB 310 
370 IF N<4 THEN L=3 @ GOSUB 310 


JPC S7 Page 28 


380 
390 
400 
410 
420 
430 
440 
450 
460 
470 


IF N<é THEN 

IF N<8 THEN 

IF N<12 THEN 
IF N<14 THEN 
IF N<18 THEN 
IF N<20 THEN 
IF N<24 THEN 
IF N<30 THEN 
IF N<32 THEN 
L=37 à GOSUB 


L=5 à GOSUB 310 
L=7 à GOSUB 310 
L=11 à GOSUB 310 
L=13 à GOSUB 310 
L=17 à GOSUB 310 
L=19 à GOSUB 310 
L=23 à GOSUB 310 
L=29 à GOSUB 310 
L=31 à GOSUB 310 
310 à GOTO 120 


eh eo ee eee eee eee ee ee eee eee ee RH HO HR RAR AR 


Programme "TRI" (tri de fichiers texte, nécessite EDLEX et JPC Rom) 


10 INPUT ‘Fichier ? !';F$ 
20 F=FILESZR(F$) 
30 ASSIGN #1 TO F$ 


40 


p=1 


50 WHILE P<IP((F+1)/9) à P=3*P+1 @ END WHILE 


60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 


REPEAT 
FOR K=0 
FOR 


NEXT 
NEXT K 
P=IP(P/3 
UNTIL P=0 
T=TIME-T 


TO P-1 
I=K+P TO F-1 STEP P 
READ #1,1;S$ 
J=I-P : 
LOOP 
READ #1,J:T$ 
IF J<K+P OR S$>=T$ THEN LEAVE 
REPLACE #1,J+P;T$ 
J=J-P 
END LOOP 
READ #1,K;U$ 
IF U$>S$ THEN J=K-P à REPLACE #1,K+P;U$ 
REPLACE #1,J+P;S$ 
I 


) 


250 ASSIGN #1 TO * 


2 3e 2 AB eee ee ee ee he ee he ee ee re eee he ee eo ee ee eee ee he he He HR HA ROAD RCRR 


Programme "TXTPC71" (conversion de fichiers en provenance de PC, Nécessite JPC Rom) 


30 


40 
50 
60 


TXTPC71 conv 
Ne traite pa 


ersion caracteres speciaux fichiers TXT PC<>TXT 71 
s Les ESC. (c) EG 06/04/88 V1.0 


DIM R1$192] ,R2$ [92] ,R3$ [96] ,R4$ [96] 


ISTART': INPUT 'Xform text: !;F1$ @ F1$=UPRC$(F1$) 


R9=FILESZR(F 
DISP 'PC>HPT 


1$) à IF R9<=0 THEN BEEP à GOTO 40 
1] or HP>PC[2]?! 
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80 
90 
100 
110 
120 
130 


140 


150 


160 


170 
180 
190 
200 
210 
220 
230 


240 
250 
260 
270 
280 
290 
300 


310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
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KS=UPRC$(KEYWAIT$) 
IF NOT POS(!12!,K$) THEN BEEP à GOTO 70 
K=VAL(K$) @ ASSIGN #1 TO F1$ 

GOSUB 'TABLES'! 

ON K GOSUB "PC2HP!,'HP2PC! 

BEEP 400,.5 à BEEP 300,.5 à DISP F1$&' done! 
ASSIGN #1 TO "!! @ GOTO START! 


IPC2HP': MAP #1,R1$,R2$ @ RETURN 
'HP2PC': MAP #1,R3$,R4$ @ RETURN 
"TABLES! : 


1T11!: ! Conversion PC>HP 
RIS=HTAS(1518191A808182838485868788898A8B8C8D8E8F ! , 1) 
R2$=HTAS('BD5E563EB4CFC5COCCC8D4B5C1CDCSDDD1D9D8D0! ,1) 
RIS=RISRHTAS( 909192939495969798999A9B9CID9IE9FAUAÏAZ! ,1) 
R2S=R2$8HTAS( 'DCD7D3C2CECAC3CBEFDADBBFBBBC20BEC4D5C6! ,1) 
RIS=RISEHTAS( !'AZA4ASAGATABAGAAABACADAEAFE1F1!,1) 
R2S=R2SRHTAS( 'C7B7B6F9FAB92020F8F7B8FBFDDEFE! ,1) 


11121: ! Equiv. appr. semi-graphiques PC 

RIS=RISSHTAS( 'BOB1B2B3B4B5B6B7B8B9BABBBCBDBERF ! ,1) 
R2$S=R2$S8HTAS('20FCFC7C2B2B2B2B2B2B2B2B2B2B2B2B! ,1) 
RIS=RIS&HTAS( 'COC1C2C3C4C5CÉC7CBCICACBCCCDCECF! ,1) 
R2$=R2$SHTAS( '2B2B2B2B2D2B2B2B2B2B2B2B2B3D232B! ,1) 
RI$S=RISSHTAS( 'DOD1D2D3D4D5D6D7DBEDIDADBDCDDDEDF ! ,1) 
R2S=R2F&HTAS( !2B2B2B2B2B2B2B2B2B2B2BFCFCFCFCFC! ,1) 


‘T21!: ! Conversion HP>PC 
R3$=HTAS(!AOATA2A3A4ASAGATABAGAAABACADAEAF! ,1) 
R4S=HTAS(12041414545454949202020207E55559C! ,1) 
R3S=R3SEHTAS( 'B0B1B2B3B4B5B6B7B8B9BABBBCBDBEBF! ,1) 
R4S=R4SEHTAS( ! 202020A78087ASA4ADAB209C9D159F9B! ,1) 
R3$=R3SEHTAS( !COC1C2C3C4C5CÉC7CECICACBCCCDCECF! ,1) 
R4S=RGSEHTAS( !83889396A082A2A3858A959784899481! ,1) 
R3$=R3SBHTAS(!FOF1F2F3F4FSFÉF7F8FOFAFBFCFDFEFF! ,1) 
R4S=RASEHTAS( ! BFBC4F9286A16F918E69999A908BE 14F! ,1) 
R3$S=R3SEHTAS('E0E1E2E3E4ESE6E7E8E9EAEBECEDEEEF! ,1) 
R4S=RASEHTAS( 1414141444449494F4F4F6F5373555979! ,1) 
R3$=R3S&HTAS(!FOFIF2F3F4FSF6F7F8FOFAFBFCFDFEFF! ,1) 
R4S=R4SEHTAS( 2020202020202DACABAGA7AEBIAFF 120! ,1) 
RETURN 








LE COIN DES LHEX 


Comme de coutume, cette rubrique contient la liste des codes hexadécimaux des fichiers Lex 
parus ce mois-ci. 


Rappelons ce qu’est un fichier Lex : c’est un programme pour le HP-71, en assembleur, qui 
apporte de nouvelles fonctions. Celles-ci sont utilisables directement, ou dans des programmes 


Basic. 


Pour bénéficier de ces nouvelles fonctions, vous n’avez pas besoin de programmer vous-même en 
assembleur, ni de posséder un module Forth/Assembleur. 


Il suffit de recopier le petit programme basic "MAKELEX" ci-dessous, de le lancer et de recopier 
les codes du fichier Lex désiré. Quand vous avez fini, les nouvelles fonctions sont accessibles, 
après avoir éteint et rallumé votre HP-71. 


Si l'erreur "Erreur de somme" apparaît, vérifiez la ligne que vous avez introduite. 


Vous trouverez donc ci-après les Lhex de la rubrique assembleur de ce mois-ci. 


LIN LIN XFN 94002 INVLIN XFN 94003 
HORNER HORNER XFN 93010 
LEXTRACE TRACE XFN 92101 


10 CALL MLEX @ SUB MLEX Q SFLAG -1 @ PURGE AH à INPUT "Nb. d'octets: ";N à LC OFF 

20 CREATE DATA AH,1,N-4 @ A=HTD(ADDR$('"AH")) à B=A à GOSUB 130 

30 Q=1 à X=0 à INPUT "000: ",P$;A$ à C$=A$ à S=0 à GOSUB 90 

40 Q=2 à X=1 à GOSUB 80 à AS=AS$&CS à A=A+37 Q N=N*2+37 à Q=3 Q SFLAG 5 @ FOR X=2 TO N DIV 16-1 
50 GOSUB 80 à C$=C$[IS*FLAG(5)+1] à POKE DTH$(A),C$ @ A=A+16-5*FLAG(5,0) @ NEXT X à Q=4 
60 DISP DTH$(X)[31; à INPUT ": 1,P$[1,MOD(N,16)]1;C$ à GOSUB 90 

70 POKE DTH$(A),C$ @ POKE DTH$(B),A$ à CFLAG -1 à END 

80 DISP DTH$(X) [3]; à INPUT ": ",P$;C$ 

90 DISP DTH$(X) [3]; @ INPUT ! sm !,11---1:D$ 

100 M=S à FOR Z=1 TO LENC(C$) à M=NUM(C$[Z]1)+M+1 @ NEXT Z 

110 IF D$=DTH$(MOD(M,4096)) [3] THEN GOSUB 130 à S=M à RETURN 

120 DISP "Erreur de somme" à BEEP à P$=C$ à POP à ON Q GOTO 30,40,50,60 

AOL PS= Se se menmnenn er " @ RETURN 
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LIN 229 octets LEXTRACE 221 octets 


HORNER 209 octets 
0123456789ABCDEF sm 0123456789ABCDEF sm 
0123456789ABCDEF sm 

000: C494E40202020202 357 000: c454854525143454 362 
001: 802E000000000000 686 000: 84F425E454250202 36A 001: 802E000000000000 691 
002: EC100E5203000000 9DE 001: 802E000000000000 699 002: FB100C5565600000 9F8 
003: F020000000000000 D06 002: 6A100D5A0A000000 9FC 003: F710000000000000 D26 
004: 013000FF0016000F 063 003: F710000000000000 D2A 004: 0D1000F945251434 086 
005: B94E465C494E4305 3FB 004: OC1000FB84F425E4 OBC 005: 54561FF621141113 3F8 
006: C494E4201FF88888 79F 005: 5425A01FF8A228FD 469 006: 7109137AF2AF014B 78B 
007: 5570D07CB01377C5 B2A 006: F8E08FFB6CO8F9D3 849 007: 2030E9029230cC906 AFO 
008: 011BAF711A8FCA4C EE8 007: D00417FD0E42030B BD3 008: 7D7BB0104A1471377 E76 
009: 077807E80799017F 269 008: A851577901118059 F41 009: EDO8F499C08D612F 231 
00A: 6D11888885577901 5DF 009: 01908D91FB00A816 2C9 00A: 07C90c6104147137 5A6 
008: 7F137742011BAF71 967 00A: 81620B0AD880DF17 663 00B: 7DB08F499C004108 937 
00C: 1A8FA34C07F40765 DO5 00B: F1331FO0C8F2141AF A07 00C: 11913717A1471371 C8E 
00D: 071607D5064E0135 06B 00C: OF6DA81681681681 D8E 00D: 7F11A81ECECE10B7 056 
O0E: 8FDFB8EO0AFB135AF6 44F 00D: 6F610BD7AF2DB892 150 00E: E808F499CO04AF71 3FD 
00F: 8FDF8E08F439C080 80D 00E: BOAFOB64540C0C98 4F5 O0F: 19AFF109DBDAD220 7C1 
010: DF89080206500298 B80 O0F: FBBCE012B8168168 8A2 010: 31218F7C210DF135 B46 
011: OCF8F363C001102A FOB 010: 168168161351C887 COB 011: 5BODB8DA93901101 ED8 
012: F910B011F0C8F214 29D 011: 0501C67E7012BA7E F9F 012: 51710CF11115171C1 245 
013: 31310117F17F011C 607 012: A7E51115378FFB6C 352 013: 31E014D8218DC32F 5DC 
014: F1CF01726079DF17 9A9 013: 0684012B058FFB6C 6ED 014: 017115B417415F38 944 
015: F137AF713778DF13 D50 014: O8F5F3D08FA34C0A AA7 015: 148A651D0810E210 CB3 
016: 7787F73BF7ABF75C 124 015: F9AF7AF67C308FFB E94 016: B173E6E6018D6CC7 066 
017: F137AF713778BF13 4C9 016: 6C08F363C012B04A 224 017: 11537AF1058FFB6C 413 
018: 7785F11BAF711A8F 870 017: 7E0512B59C1FO0C8F 5D3 018: 00411BCE4A410B7E 7AB 
019: CA4CO738F7A8F137 C31 018: 21471358218F499C 94C 019: 201577AF3058F049 B2B 
O1A: AF71377C8F011371 FBF 019: 08D612F0870B017F CD6é O1A: CcO8F363C0731011B EAA 
01B: FOC8F2145137018F 34C 01A: 15370117828AF015 03E 01B: CE10B5CD11913701 238 
01C: DF8E08FFB6CO11BA 727 01B: 3181481481481481 38F 01C: 11A04F2D7137CB13 5C6 
01D: F711A8F363C08F49 ACF 01C: 481401F 50E 01D: 701048228216CDF 90B 


O1E: 9C0704F8F612F0F E3E 
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ENGLISH SUMMARY 
JPC 57 - SEPTEMBER 1988 


We are very happy to report that JPC Rom is now distributed by CMT worldwide. We have been very proud 
to see the advertisement in the latest EduCALC catalog. 


We are sorry to report to our foreign friends that, from now on, we will no more be able to accept 
Eurochèques payments. This is due to very high fees taken by our french banks. So, you can pay your 
subscription by the following means exclusively : 

- Postal Transfer to our account 18 823 40 C in Paris, or 

- checks in French Francs drawn on a French bank. 


The HP-28 columns begins with an article by Philippe Heïlbronn on curve fitting. It provides the best fit 
among the classical four types. Sébastien Laglande gives us an example of HP-28 assembly language 
programming with a little program to exercise the beeper. Next, Paul Courbis and Sébastien explain the 
solution of their game published in JPC 54. Last, Paul gives a shorter version of the INKEY program published 
in JPC 56. 


The first HP-75 article is from Eric Gengoux : it explains the various possibilities to exchange data files 
between HP-75 VisiCalc and PC compatible using V.C.S. format. Then, Robert Atlan gives us a new 
method to compute linear regressions using Tukey method. 


The HP-71 assembly language columns begins with a Lex file from Jean Dardennes to help you compute 
linear interpolation. Next, Laurent Istria gives an answer to Jack Elhay (cf JPC 50) with a new version of the 
Horner method to compute the value of a polynom. Last article in this columns is from Lionel Guillou and 
provides a keyword to compute the sum of the diagonal elements of a square matrix. 


A new basic version of the well known prime number determination problem by Philippe Aspero. Then, 
Laurent Istria studies various sorting methods of text files, and last, Eric Gengoux provides us an example of 
use for MAPS function to convert PC compatible characters to Roman 8. 


Until next month, 


Happy Programming and JPC reading ! 


